AD/A-001  387 


DYNAMIC  ANALYSIS  OF  A  BUILDING  AND  BUILDING  ELEMENTS 


Stanford  Research  Institute 


PREPARED  FOR 

Defense  Civil  Preparedness  Agency 
April  197^ 


DISTRIBUTED  BY: 


National  Toclmical  lafanMtion  Sirvico 
U.  S.  DEPARTMENT  DF  COMMERCE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Whan  Oat*  Entarad) 


'  AfAO  INSTRUCTIONS 
SSFOne  COMPLCTING  FORM 


3.  OOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

(none) 


4.  TITLE  land  SLibtltla) 


DYNAMIC  ANALYSIS  OF  A  BUILDING  AND 
BUILDING  ELEMENTS 


8.  TYPE  CF  REPORT  S  PERIOD  COVERED 

Technical  Report 


8.  PERFORMING  ORO.  REPORT  NUMBER 


7.  AUTHORItl 


C.  K.  Wlehle 


B.  CONTRACT  OR  GRANT  NUMBERIol 

DAHC20-71-C-0292 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 
Stanford  Research  Institute 
Facilities  and  Housing  Research  Dept. 
Menlo  Park,  California  94025 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Defense  Civil  Preparedness  Agency 
Washington,  D.C.  20301 


14  MONITORING  AGENCY  NAME  ft  ADDRESS  (If  dlH.  from  Controlling  Offical 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

DCPA  Work  Unit  11541 


12.  REPORT  DATE  13.  NO.  OF  PAGES 

April  1974  /5/ 


18.  SECURITY  CLASS,  (of  thit  raport) 


UNCLASSIFIED 


16.  DISTRIBUTION  STATEMENT  (of  tfiit  raport) 


ISa.  DECLASSIFICATION /DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  th*  abttraet  antarad  in  Block  30,  if  diffarant  from  raport) 


IB.  SUPPLEMENTARY  NOTES 


WICES  SUSJEa  TO  CHANGE 


19.  KEY  WORDS  (Continu*  on  ravatM  tkf*  if  nacaiaary  and  Idantify  by  block  number) 

Restrained  reinforced  concrete  slabs 
Wood-Joist  floor  systems 

Dyn^lc  analyst,  "•"NATIONAL  TECHNICAL 

INFORMATION  SERVICE 

r  rames  JJ  g  Depailment  of  Commerce 

VA  ZZ151 


30.  ABSTRACT  I  Continue  on  ravari*  (Id*  If  nacataary  and  klcnr'fy  by  block  number) 

The  objective  of  the  overall  rebearch  program  is  to  develop  an  evaluation 
procedure  applicable  to  existing  NSS-type  structures  and  private  homes.  Past 
efforts  have  been  concerned  with  examining  exterior  walls,  window  glass,  steel 
frame  connections,  applications  to  actual  buildings,  and  reinforced  concrete 
fljor  systems.  This  report  presents  modifications  to  the  mathematical  model 
previously  developed  fur  reinforced  concrete  floor  slabs  to  include  the  effect 
of  longitudinal  restraint  on  slab  resistance.  Also  included  in  the  report  are 

V 


DD.r'T.wa 

EDITION  OF  1  NOV  M  II  OMOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLABSIPICATION  OF  THIS  PAGE  (Whan  Data  Entarad) 


«CumTY  CLAMIFICATIQN  Of  TMII  PAQt  (Wl«n  QlM  tnwn) 

19.  KEY  WORDS  (Contlnutd) 


30  abstract  (Continu«dl 

the  development  of  the  resistance  function  for  predicting  the  collapse  of  wood- 
joist  floor  systems,  and  the  dynamic  inelastic  analysis  of  an  eight-story  steel- 
frame  office  building. 


DD.  XaMTa*""' 

EDITION  Of  1  NOV  ES  IS  OBSOLETE 


iii _ 

_ UNCLASSIFIED 


SECURITY  CLASSIf  ICATION  OF  THIS  PAGE  (Whan  0*M  Entarad) 


STANFORD  RESEARCH  INSTITUTE 

Menlo  Park  California  9402S  USA 


I 


Technical  ftaport 


April  1974 


DYNAMIC  ANALYSIS  OF  A  BUILDING 
AND  BUILDING  ELEMENTS 


By:  C.  K.  WIEHLE 

FaeiUtitt  and  Homing  Hataareh 


Prapared  for: 

DEFENSE  CIVIL  PREPAREDNESS  AGENCY 
WASHINGTON,  D.C.  20301 


CONTRACT  DAHC20-71-C-0292 
DCPA  Work  Unit  11541 


SRI  Projea  1219-2 


Approved  for  public  release;  diitrlbution  unlimited. 


DCPA  REVIEW  NOTICE 

Thii  report  hat  been  reviewed  in  the  Defeneo  Civil  Preparedneii  Agency  and  approved  for  publication. 
Approval  doei  not  tignify  that  the  content!  nacanarily  reflect  the  rdViewi  and  pollciae  of  the  Dafenie 
Civil  Preparadnatt  Agency. 


•  I 

II 


ABSTRACT 


The  objective  of  the  overall  research  program  is  to  develop  an 
evaluation  procedure  applicable  to  existing  NSS-type  structures  and 
private  homes.  Past  efforts  have  been  concerned  with  examining  exterior 
walls,  window  glass,  steel  frame  connections,  applications  to  actual 
buildings,  and  reinforced  concrete  floor  systems.  This  report  presents 
modifications  to  the  mathematical  model  previously  developed  for  reinforced 
concrete  floor  slabs  to  include  the  effect  of  longitudinal  restraint  on 
slab  resistance.  Also  included  in  the  report  are  the  development  of  the 
resistance  function  for  predicting  the  collapse  of  wood-Jolst  floor 
systems,  and  the  dynamic  inelastic  analysis  of  an  eight-story  steel-frame 
office  building. 
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SUMMARY 


Introduction 

Tho  objoctlvu  of  thii  Invest  igntloa  was  to  develop  nn  evslustion 
procedure  for  d«teraining  the  blast  protection  afforded  by  existing  NSS- 
type  structures  and  private  residences.  The  approach  adopted  was  to 
fonailate  a  procedure  for  exaaining  the  dynaaic  response  of  a  structure 
over  a  range  of  incident  overpressure  levels  to  detemine  the  jtpsseiins 
at  which  ooll^MO  of  the  various  elenents  occurs. 

Background 

Past  efforts  in  this  progran  have  exanined  exterior  walls,  window 

.4 

glass,  steel-fraae  connections,  applications  of  the  procedures  to  actual 
NSS  buildings,  and  reinforced  concrete  floor  systens.  Hie  report  sun- 
narises  the  three  phases  of  the  research  effort  conducted  during  the 
eurront  contract  period.  First,  are  the  nodiflcations  to  the  previously 
developed  analytical  expressions  for  reinforced  concrete  floor  systens  to 
include  the  effect  of  longitudinal  edge  restraint  on  the  resistance  of  the 
slab.  (An  appendix  presents  the  correlations  between  the  analytical  pre¬ 
dictions  and  published  experinental  data  on  the  response  of  longitudinally 
restrained  slabs.)  Second,  the  report  contains  the  developnent  of  the 
resistance  function  used  to  predict  the  dsmanlc  response  and  collapse  of 
wood- joist  floor  systeaw.  Third,  a  number  of  elastic  and  inelastic  analyses 
were  performed  to  determine  the  response  of  an  elght-st(»ry,  steel-frsM 
office  building  subjected  to  nuclear  air  blast.  The  results  of  the  frame 
analyses,  together  with  an  estimate  of  the  collapse  strength  of  the  build¬ 
ing,  are  presented  in  the  report. 


Discassion 


Longitudinally  lte»tr>liM>d  Reinforced  Concr»te  Slaba 

Since  there  was  a  nded  to  determine  the  blast  strength  of  floor 
systems  for  actual  MFSS  buildings,  the  previous  effmrt  was  concerned 
primarily  with  the  development  of  Interim  analytical  procedures  for  pre¬ 
dicting  the  dynamic  response  and  collapse  of  various  types  of  reinforced 
concrete  floor  systems  over  basement  areas.  During  the  current  study. 

It  was  possible  to  modify  the  resistance  function  for  floors  by  Including 
the  effect  of  longitudinal  edge  restraint  on  the  slab  resistance.  Com¬ 
pressive  membrane  forces,  which  occur  In  slabs  that  are  restrained  from 
outward  movement  at  their  edges,  can  be  a  controlling  factor  In  determining 
the  magnitude  of  the  slab  resistance;  therefore.  Its  Inclusion  In  the 
mathematical  model  enhances  the  ability  to  make  realistic  collapse  pre¬ 
dictions  for  blast-loaded  flocrs  In  existing  buildings.  It  was  also 
possible  to  make  comparisons  between  analytical  predictions  and  various 
published  experimental  data  on  the  static  and  dynamic  response  of  re¬ 
strained  slabs.  The  results  of  the  analytical  and  experimental  correlations, 
presented  In  the  report.  Justify  the  use  of  the  compressive  membrane  mode 
for  predicting  the  collapse  of  floor  systems  in  existing  structures  that 
are  longitudinally  restrained  at  the  edges. 

Wood-Joist  Floor  Systems 

One  of  the  goals  of  this  research  program  has  been  to  examine  the 
blast  strength  of  private  residences.  In  previous  studies,  analytical 
procedures  were  developed  for  predicting  the  collapse  of  exterior  walls  of 
typical  residences,  and  the  analytical  predlctl<aiB  were  verified  by  com¬ 
parison  with  results  of  nuclear  weapon  field  experiments.  Because 
of  the  possibility  of  using  home  basements  for  blast  shelters  during 
emergencies,  the  development  of  a  method  to  predict  the  collapse  of 


wood-Jolst  floors  was  initiated  during  the  current  research  study.  A 
literature  search  of  the  pertinent  publications  Indicated  that  wood-Jolst 
floor  systens  are  complex  structural  systems  that  are  not  readily  amenable 
to  precise  theoretical  analysis  because  of  the  wide  variation  In  the 
mechanical  properties  of  wood,  and  the  variability  In  construction  methods 
and  details. 

The  resistance  function  developed  for  wood-Jolst  floors  consists  of 
elastic  and  Inelastic  phases.  The  elastic  phase  Is  limited  by  the  bending 
strength  of  the  weakest  Joist  (first  break),  which  Is  limited  by  a  defect 
resulting  from  the  largest  permissible  knot  for  the  grade  and  species  of 
wood  used  In  the  Joists.  The  Inelastic  phase,  which  results  from  the 
nonlinear  behavior  of  wood  In  compression  above  the  proportional  limit. 

Is  limited  by  the  theoretical  ultimate  strength  of  a  clear  wood  specimen. 
The  failure  criterion  Is  developed  In  the  form  of  a  truncated  normal 
probability  function. 

Dynamic  Analysis  of  North  Carolina  National  Bank  Building 

A  continuing  concern  In  evaluating  the  collapse  overpressure  of 
existing  buildings  has  been  the  relative  blast  strength  of  the  exterior 
walls  and  frames  of  multistory  buildings.  It  Is  often  assumed  for  the 
analysis  of  blast-loaded  frame  buildings  that  the  exterior  walls  can  be 
oonsldered  as  frangible,  and  therefore,  that  the  wall  loading  transferred 
to  the  frame  can  be  approximated  by  an  Impulse  loading.  However,  for 
many  of  the  actual  buildings  analyzed,  the  strength  nt  the  exterior  walls 
Wider  blast  loading  was  sufficiently  high  to  make  It  doubtful  that  the 
frame  could  survive  at  the  overpressure  level  required  to  collapse  the 
walls.  Therefore,  In  this  phase  of  existing  oomputer  program  for  analy¬ 
sing  the  elastic  and  Inelastic  dynamic  response  of  two-dimensional  struc¬ 
tural  frames  was  used  to  analyse  the  frame  of  an  eight-story  steel-frame 
bank  and  office  building.  This  building  was  selected  for  study  because 


a  previous  analysis  had  Indicated  that  the  exterior  walls  had  an  Incident 
collapse  overpressure  of  15.7  psl,  and  the  building  could  therefore  be 
considered  as  strong-walled. 

A  series  of  elastic  and  plastic  analyses  were  performed  to  determine 
the  dynamic  response  of  the  frame  of  the  building  to  various  blast  over¬ 
pressure  levels.  Even  though  the  computer  program  used  could  not  predict 
frame  collapse,  the  results  of  the  analyses  Indicated  that  the  frame 
would  probably  collapse  between  3-  and  4-psl  Incident  overpressure  level. 
For  the  case  examined,  the  blast  resistance  of  the  building  frame  was 
predicted  to  be  much  less  than  that  of  the  exterior  walls,  which  Is  an 
Important  consideration  In  predicting  either  the  magnitude  of  building 
damage  or  the  number  of  casualties  that  might  occur  in  buildings  subjected 
to  nuclear  air  blast. 


I  INTRWUCriON 


Under  contract  to  the  Defense  Civil  Preparedness  Agency,  Stanford 
Research  Institute  is  developing  a  procedure  for  the  evaluation  of  existing 
structures  subjected  to  nuclear  air  blast.  The  objective  of  the  overall 
program  Is  to  develop  an  evaluation  procedure  applicable  to  existing  NSS- 
type  structures  and  private  homes.  The  purpose  of  this  phase  of  the  re¬ 
search  was  to  improve  the  prediction  techniques  for  reinforced  concrete 
floor  systems,  develop  analytical  procedures  for  dynamically-loaded  wood- 
joist  floor  systems,  and  make  a  preliminary  examination  of  the  response  of 
multistory  buildings  to  nuclear  blast  loading. 

Background 

Past  efforts  in  this  program  have  been  concerned  with  examining 
exterior  walls  (Ref.  1),  window  glass  (Ref.  2),  steel-frame  connections 
(Ref.  3),  two-way  action  walls  (Ref.  4),  applications  to  NSS  buildings 
(Refs.  5  and  6),  and  reinforced  concrete  floor  systems  (Ref.  7). 

Even  though  the  primary  Interest  since  the  inception  of  this  program 
has  been  in  the  behavior  and  collapse  of  the  building  system,  the  com¬ 
plexity  of  an  overall  building  evaluation  procedure  necessitated  the 
establishment  of  a  ranking  of  structural  elements.  It  is  apparent  that 
the  collapse  of  the  exterior  walls  of  most  buildings  is  important  to  the 
casualties  produced.  This  is  especially  true  for  large  multistory  build¬ 
ings  where  the  collapse  of  the  exterior  and  interior  walls  could  result 
in  a  large  number  of  casualties  as  a  result  of  ejection  from  the  building, 
even  if  the  floors  and  frame  remained  intact.  Also,  for  a  load-bearing- 
wall  building,  where  the  exterior  walls  are  the  primary  structural  member. 
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the  collapse  of  the  exterior  walls  can  precipitate  a  catastrophic  col¬ 
lapse  of  the  entire  building.  Since  one  of  the  primary  uses  of  a  build¬ 
ing  evaluation  procedure  Is  to  provide  Input  for  predicting  the  survival 
of  people  located  In  buildings  subjected  tc  nuclear  air  blast,  the 
Initial  research  effort  was  directed  toward  development  of  a  method  for 
detenpinlng  the  response  and  collapse  of  exterior  wall  elements.  The 
three  basic  types  of  exterior  walls  considered  were  unreinforced  concrete 
or  masonry  unit  walls  without  arching,  unreinforced  concrete  or  masonry 
unit  walls  with  arching,  and  reinforced  concrete  walls.  With  minor 
modification,  the  analytical  procedures  for  exterior  walls  also  apply 
to  Interior  partitions  of  similar  construction. 

As  noted  in  Ref.  7,  there  was  a  need  to  develop  Interim  procedures 
for  predicting  the  collapse  of  floors  over  basement  areas  In  existing 
NSS  buildings.  Therefore,  mathematical  models  for  the  behavior  of 
reinforced  concrete  floor  systems  were  formulated,  based  on  readily 
available  technical  Information.  During  the  current  research  effort  It 
was  possible  to  modify  the  previously  developed  analytical  expressions 
to  Include  the  effect  of  longitudinal  edge  restraint  on  the  resistance 
of  reinforced  concrete  slabs.  Compressive  membrane  forces  occur  In  con¬ 
crete  slabs  that  are  restrained  from  outward  movement  at  the  edges  as  a 
result  of  adjacent  floor  panels,  or  such  other  restraints  as  heavy 
spandrel  beams.  Since  compressive  membrane  forces  can  be  a  controlling 
factor  In  determining  the  magnitude  of  slab  resistance.  Its  inclusion  In 
the  floor  prediction  schemes  enhances  the  ability  to  make  realistic 
collapse  predictions  for  floors  In  existing  buildings. 

Although  collapse  predictions  have  been  made  for  blast-loaded  wood- 
joist  floors  on  an  Individual  basis  In  this  program,  the  procedure  was 
not  systematized.  During  the  current  effort,  a  large  body  of  analytical 
and  experimental  Information  on  wood  and  wood-joist  floors  was  reviewed 
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and  a  resistance  function  developed  for  modelling  the  complex  structural 
behavior  of  wood-joist  floors. 

A  continuing  concern  in  this  program  has  been  the  relative  strength 
of  the  exterior  walls  and  frames  of  multistory  buildings.  To  investigate 
the  relative  strength  of  the  walls  and  frame  would  require  a  comprehensive 
computer  program  that  included  the  inelastic  response  of  the  frame  under 
dynamic  loading^  as  well  as  realistic  frame  collapse  mechanisms.  Since 
no  such  program  could  bo  developed  during  this  study,  an  available 
computer  program  for  analyzing  the  elastic  and  inelastic  dynamic  response 
of  two-dimensional  structural  frames  (Ref.  8)  was  used  to  estimate  the 
blast  strength  of  an  eight-story  steel-frame  office  building.  Although 
the  frame  program  does  not  include  frame  collapse  mechanisms,  the  results 
of  the  analyses  provided  considerable  insight  into  the  relative  collapse 
strengths  of  the  exterior  walls  and  frame  of  a  typical  multistory 
building. 

The  analytical  method  used  in  the  research  study  has  been  to  establish 
the  resistance  function  for  each  structural  element  of  interest  by  con¬ 
sidering  the  approximate  response  mode,  and  by  assuming  that  the  element 
is  subjected  to  a  uniformly  distributed  static  load.  To  analyze  the 
dynamic  response  and  collapse  of  the  member,  it  is  transformed  into  an 
equivalent  single-degree-of-freedom  system  by  the  use  of  transformation 
factors  for  the  load,  resistance,  and  mass.  The  equation  of  motion  is 
then  solved  on  a  computer  using  a  numerical  integration  procedure.  Al¬ 
though  established  analytical  procedures  have  been  used  wherever  possible, 
it  has  been  necessary  to  modify  and  adapt  current  procedures  as  well  as 
develop  new  procedures  for  specific  use.  Relatively  simplified  analytical 
models  have  been  used  for  wall  and  floor  element  analysis  to  prevent  the 
evaluation  of  a  structure  from  becoi  unwieldly  as  a  result  of  excessive 
computation. 
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As  discussed  In  Ref.  4,  the  analysis  of  actual  building  elements 
subjected  to  nuclear  air  blast  requires  the  assumption  of  values  for  many 
of  the  physical  properties  of  the  structure  that  are  unknown  and  cannot 
be  measured  without  an  unwarranted  amount  of  effort.  Similarly,  assump¬ 
tions  are  also  required  in  the  determination  of  the  parameters  defining 
the  load  acting  on  the  building  element.  Since  precise  values  cannot 
usually  be  specified  for  many  of  the  parameters  that  influence  the  collapse 
of  actual  structures,  a  probabilistic  approach  was  formulated  to  provide 
a  more  realistic  method  for  evaluating  existing  structures  subjected  to 
air  blast  loading. 

Report  Organization 

Section  II  gives  the  mathematical  expressions  developed  for  modifying 
the  resistemce  functions  for  reinforced  concrete  slabs  to  Include  the 
effect  of  longitudinal  edge  restraint.  The  correlation  of  the  analyti "al 
predictions  with  experimental  data  on  the  response  of  longitudinally 
restrained  slabs  is  Included  in  the  Appendix.  The  development  of  the 
resistance  function  for  predicting  the  behavior  of  wood-joist  floor  systems 
is  presented  in  Section  III.  The  results  of  the  elastic  and  Inelastic 
analyses  of  the  dynamic  response  of  the  multistory  office  building  is 
given  in  Section  IV. 
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II  LONGirUDINALLY  RESTRAINED  REINFORCED  CONCRETE  SLABS 

Introduction 

Interim  analytical  procedures  for  calculating  the  resistance  function 
for  reinforced  concrete  slabs  Is  presented  In  Ref.  7.  In  that  study; 
the  yield-line  theory  was  used  to  predict  the  ultimate  flexural  strength 
of  the  slab.  For  slabs  with  continuous  reinforcement;  a  tensile  membrane 
mode  was  also  Included  in  the  resistance  function;  but  such  refinements 
as  compressive  membrane  action  In  slabs  with  restrained  edges  were  not 
Included  In  the  interim  techniques. 

However,  several  Investigators  (c.g..  Refs.  9,  10;  11;  and  12)  have 
found  that  the  yield-line  theory  grossly  underestimates  the  strength  of 
slabs  with  longitudinal  edge  restraint;  such  as  that  provided  by  surround¬ 
ing  floor  panels,  or  by  heavy  spandrel  beams  or  walls.  When  a  slab  with 
edge  restraint  deflects  under  load,  significant  compressive  membrane 
forces  are  Induced  In  the  slab  that  result  In  ultimate  moments  at  the 
yield  lines  considerably  higher  than  when  the  forces  are  not  present. 

Since  the  ultimate  strength  predicted  by  the  yield-line  theory  Is  too 
conservative  for  longitudinally  restrained  slabs,  the  resistance  function 
for  reinforced  concrete  slabs  presented  In  Ref.  7  has  been  modified  to 
Include  the  compressive  membrane  mode.  The  analytical  procedure  developed 
in  this  report  for  calculating  the  resistance  function  of  longitudinally 
restrained  slabs  was  based  primarily  on  the  Investigations  presented  In 
Refs.  9  and  10. 

As  shown  on  Figure  1,  the  resistance  function  for  longitudinally 
restrained  slabs  consists  of  three  phases: 
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RESISTANCE 


FIGURE  1  RESISTANCE  FUNCTION  FOR  LONGITUDINALLY 
RESTRAINED  REINFORCED  CONCRETE  SLAB 
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•  Compressive  membrane  resistance 

•  Secondary  resistance 

•  Tensile  membrane  resistance. 

The  development  of  the  analytical  expressions  for  determining  the  slab 
resistance  for  each  of  these  phases  Is  presented  In  the  following  sub¬ 
sections.  The  analytical  methods  for  predicting  the  ultimate 

flexural  strength  for  the  compressive  membrane  resistance  phase  are  based 
primarily  on  the  Investigations  summarized  In  Refs.  9  and  10.  The  methods 
used  for  determining  the  slab  resistance  during  the  secondary  and  tensile 
membrane  resistance  phases  are  based  on  the  expressions  previously  de¬ 
veloped  in  this  program  (Ref.  7). 


Compressive  Membrane  Resistance 


The  ultimate  resistance  during  the  compressive  membrane  phase  Is 
determined  by  assuming  the  slab  to  fall  along  the  yield  lines  shown  In 
Figure  2.  The  value  of  defining  the  location  of  the  yield  lines,  Is 

determined  from  the  fallowing  equation,  previously  given  in  Ref.  7: 


1  \g  Va 

2  vi 


(1) 


where 


^  2  Vl  + 

(la) 

=  2  Vl  ^ 

• 

(lb) 

Each  segment  of  the  pattern  Is  assumed  to  rotate  as  a  plane  surface. 

For  a  slab  with  edges  restrained  against  longitudinal  movement,  both 
thrust  and  moment  act  on  sections  along  the  yield  lines.  These  thrust 
forces  significantly  Increase  the  moment  resistance  of  the  slab  cross 
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FIGURE  2  ASSUMED  COLLAPSE  MECHANISM  FOR 

LONGITUDINALLY  RESTRAINED  REINFORCED 
CONCRETE  SLAB 
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section.  This  Increase  is  determined  by  using  the  moment  equations 
developed  previously  In  Appendix  A  of  Ref.  1.  In  this  study,  however, 
the  moments  are  related  to  the  mid-thickness  of  the  slab  rather  than  to 
the  centroid  of  the  tensile  reinforcement.  This  Is  done  since,  as  sub¬ 
sequently  shown,  the  thrust  forces  are  also  assumed  to  act  at  the  mid¬ 
thickness  of  the  slab. 

The  magnitude  of  the  thrust  induced  In  the  longitudinally  restrained 
slab  when  the  collapse  mechanism  forms  is  determined  by  assuming  the  slab 
to  be  composed  of  strips  oriented  In  the  x  and  s  directions,  as  shown  in 
Figure  3.  A  strip  In  either  direction  has  the  same  thickness  as  the  slab 
and  contains  only  the  reinforcement  steel  in  that  direction.  The  yield 
sections  are  perpendicular  to  the  direction  of  the  strip,  with  the 
torsional  moments  on  these  sections  assumed  to  be  zero.  The  static  load 
resistance  can  now  be  expressed  as  the  sum  of  the  load  carried  by  each 
of  the  individual  strips. 

Consider  a  typical  strip  4-2-2 '-4 '  in  the  x-direction,  with  plastic 
hinges  Just  formed  at  4,  2,  2^,  and  4^,  as  shown  in  Figure  4.  The  forces 
acting  on  portion  4-2  at  this  stage  of  behavior  are  shown  in  Figure  5. 

In  addition  to  the  yield-line  moments,  compressive  membrane  forces  are 
induced  by  longitudinal  restraint  at  the  ends  of  the  strip.  Considering 
the  equilibrium  of  these  horizontal  forces, 

Ka  =  =  Nu,  .  (2) 

The  subscript  u  refers  to  the  ultimate  stage  of  behavior  (initial  crushing 
of  the  extreme  fiber  of  the  concrete  at  each  hinge  section). 

The  strain  at  the  extreme  fiber  of  the  concrete  is  assumed  to  be 
equal  to  the  crushing  strain  of  the  concrete.  Assuming  the  crushing 
strain  is  reached  simultaneously  at  sections  2  and  4,  the  thrust  at 
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FIGURE  3  SYSTEM  OF  EQUIVALENT 
SLAB  STRIPS 
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FIGURE  4  COLLAPSE  MECHANISM  OF  A  TYPICAL  STRIP  IN  x-DIRECTION 


FIGURE  5  INTERNAL  FORCES  AT  YIELD  SECTIONS  OF  STRIP  4-2 


section  2  in  Figure  6  required  to  Initiate  crushing  of  the  concrete  Is 
given  by 


-  Pgdjf.e  +  ft,'daf/a 


or 


N«a  = 


where 


(3) 


(3a) 


Similarly,  the  thrust  at  section  4  (Figure  7)  required  to  Initiate  crush¬ 
ing  of  the  concrete  Is 

Nu4  =  -  ?4)  (4) 


where 


(4a) 


An  enlarged  view  of  the  strip,  relating  the  depths  to  the  neutral 
axis,  kydg  and  k„d4,  to  the  geometry  of  the  strip  at  a  deflection  y.  Is 
shown  In  Figure  8.  The  effect  of  partial  restraint  against  In-plane 
movement  at  the  supports  is  taken  Into  account  by  assuming  the  support 
lines  at  4  and  A*  to  move  outward  a  distance  s,l^/2.  This  longitudinal 
movement  Increases  the  original  distance  between  supports  from  1^  to 
(1  +  s,)I^. 

The  slab  also  shortens  axially  due  to  the  effects  of  the  compressive 
membrane  force.  The  total  axial  strain,  e^,  will  have  a  constant  value 
along  the  length  of  the  strip,  since  the  membrane  force  Is  constant. 
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FIGURE  6  CONDITIONS  AT  SECTION  2  OF  STRIP  4-2 
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FIGURE  7  CONDITIONS  AT  SECTION  4  OF  STRIP  4-2 
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Because  of  e, ,  the  shortening  of  the  middle  portion  2-2'  of  the  strip 
will  be  (l^  -  2x)  and,  because  of  symmetry,  the  ends  of  the  portion 
2-2'  will  approach  the  center  of  the  strip  by  0.5  e,  (li  -  2x).  Also,  due 
to  Sit,  the  lengths  of  portions  4-2  and  4^-2'  will  decrease  to  (1  -  e^)x. 
The  distance  between  points  4  and  2  (Figure  8)  is  given  by 

r  1 

I  X  +  —  +  0.5  €,  (li_  -  2x)  sec  cp  =  (1  -  c,  )x 

+  (hj  -  k^d4 )tan  ^ 


-  kudgtan  cp 


or 


hg  “  ^  ^  — 


nr;  [("  f)*  *  '•(*  ““  f  +  r  -  *) 

fiiil 

^  2  J  • 


Since  m  is  small  (thus  sin  ;p  ^  2  sin  ^  ^  ^  and  cos  (d^  1),  Eq.  5  reduces  to 

2  X 


y 

hg  -  =  2  “ip 


(6) 


Therefore, 


kydg  +  kydt  = 


y  xKs/ 
■>.  - 1  -  — 


(7) 


where 


Sx  =  Sx  +  Sa  .  (7a) 

From  the  geometry  of  the  deflected  shape,  shown  in  Figure  4, 
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>' '  ””  (ik) 


X  s  BLl 


(8) 


Substltu^^lng  Eq,  8  in  Eq.  7  yields 


kydg  +  kud4 


(9) 


Substituting  Eqs.  3  and  4  Into  Eq.  2  yields 


kikyda  -  Ts  kikud4  -  r4 

or 

kydg  -  kyd4  =  (r^  “  r* )  .  (10) 

•*1 


Solving  Eqs.  9  and  10  simultaneously  yields  the  following  two  expressions 
for  the  depths  to  the  neutral  axis 


kuOi 


Z  <^2  -  ^4)! 


1  _  J 

r-  <*-4  -  rg  )/ 


(11) 


(12) 


Examination  of  these  equations  shows  that  the  effect  of  s/  Is  to  reduce 
the  depths  to  the  neutral  axis  at  the  yield  sections. 


Substituting  Eq.  11  Into  Eq.  4,  and  noting  from  Eq.  2  that  =  Ny4 , 
the  thrust  developed  at  the  ends  of  the  strip  Is 


N 


MX 


(13) 


To  evaluate  from  Eq.  13,  It  Is  necessary  to  know  the  values  of  fg 
and  r^ .  These  depend  upon  the  steel  stresses  f,a ,  f/g,  f ,4  ,  and  f ,'4 . 
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The  values  of  these  stresses  In  turn  depend  on  the  value  of  the  balanced 
thrust  for  the  cross  sections  4  and  2, 


The  balanced  thrust  Is  the  thrust  that  causes  yielding  of  the  tension 
steel  and  crushing  of  the  concrete  to  occur  simultaneously.  Its  value 
Is  dependent  on  the  properties  of  the  cross  section.  For  any  thrust  less 
than  the  balanced  thrust,  N^,  the  tension  steel  Is  yielding  (f,  =  f^) 
when  crushing  of  the  concrete  occurs.  It  Is  also  assumed  that  the  com¬ 
pression  steel  Is  yielding  (f/  =  f^)  when  crushing  occurs,  as  will 
generally  be  the  case.  Referring  to  Figures  6  and  7,  the  positions  of 
the  neutral  axis  at  each  hinge  section  for  a  balanced  thrust  are 


’  [c.  t  (f./E,)]"* 

The  subscript  b  refers  to  the  balanced  condition. 


Substituting  Eqs.  14  and  15  Into  Eqs.  3  and  4,  respectively,  the 
balanced  thrusts  are  given  by 


(14) 


(15) 


»b.  =  f;>H 


(16) 


(17) 


where 


ra  =  dg  (P8 

-  Pa)  TF 

(16a) 

r*  =  <P4 

-  R4  )  ^ 

(17p) 
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For  the  case  of  a  tension  failure  at  both  ends  of  strip  4-2,  the  thrust 
at  initial  crushing  of  the  concrete  is  found  from  Eq.  13  to  be 


N 


uv 


provided 


Nu,  <  N„ 

Nu,  «  Nk4 


(18) 


(18a) 


Next,  considering  a  typical  strip  3-1-1 '-3 '  in  the  z-direction 
(see  Figures  2  and  3),  the  strips  are  similar  to  those  in  the  x-direction, 
except  that  the  strips  are  of  length  L,,  rather  than  .  Also,  across 
the  central  portion  of  the  slab  only  two  segments  exist  instead  of  the 
three  segments  across  the  diagonals.  An  enlarged  view  of  the  diagonal 
portion  of  segment  3-1  is  shown  in  Figure  9.  Comparing  this  with  Figure  8, 
the  geometry  is  seen  to  be  the  same  except  for  the  change  in  notation 
from  X,  ,  4,  and  2  for  the  x-direction  strips  to  z,  Lj ,  3,  and  1  for 
the  z-direction  strips.  The  corresponding  equations  for  the  z-direction 
strips  can  thus  be  obtained  by  merely  interchanging  the  notation.  From 
Eq.  7, 


k^di  +  k^d. 


h. 


2  "  2y 


(19) 


where 


s,'  =  s,  +  €, 


From  the  geometry  of  the  deflected  shape. 


0  S  z  S  — ■ 
2 


(19a) 


(20) 
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This  equation  is  the  sane  as  Eq.  8  if  P  =  1/2.  Substituting  this  value, 
along  with  the  previous  notation  changes,  in  Eqs.  11  and  12,  the  depths 
to  the  neutral  axis  are  given  by 


1  (r  , 

1 

-  ^3  >j 

- ;  j['‘-  - 1 

“  4y„. 

(21) 

1  (f  / 

z\  8*^11 

1 

1 

^^3  " 

■  5  tl’'-  ■  ( 

Ls)^“  "  4y„  J 

(22) 

The  values  of  r^  and  Fg  depend  on  the  value  of  the  balanced  thrust  for 
the  cross  sections  at  1  and  3.  From  Eqs.  16  and  17,  the  balanced  thrusts 
are  given  by 


where 


[e.  +  Cf,/E,)] 


d,  -  f"r 


Nbs  =  f 


da  - 


ri  -  di  (pi 

**3  =  dg  (Pa 


(23) 

(24) 

(23a) 

(24a) 


For  the  case  of  a  tension  failure  at  both  ends  of  strip  3-1,  the  thrust 
at  initial  crushing  of  the  concrete  is  found  from  Eq.  18  to  be 


(25) 


provided 


Nu,  ^  Nbi  j 

Nu*  ^  NfcS  )  .  (25a) 
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Ultimate  Deflection 


The  ultimate  deflection,  yy ,  Is  the  central  deflection  of  the  slab 
at  Initial  formation  of  the  flexural  collapse  mechanism.  In  other  words, 
Yy  Is  the  deflection  required  to  develop  the  crushing  strain,  Cy ,  along 
the  hinge  lines.  For  the  general  case  of  a  rectangular  slab,  this  value 
will  be  different  for  the  x-dlrectlon  and  z-dlrectlon  strips  due  to  the 
difference  In  the  length  of  these  strips.  Only  for  the  special  case  of 
a  square  slab  will  yy  be  the  same.  Since  the  short  direction  strips 
usually  provide  the  greatest  portion  of  the  slab  resistance,  the  approach 
taken  In  this  study  Is  to  use  the  value  of  yy  for  the  z-dlrectlon  (short 
direction)  strips,  except  for  the  cases  subsequently  noted.  This  value 
Is  calculated  as  follows. 

From  the  geometry  of  the  z-dlrectlon  strip  3-1  shown  In  Figure  9, 

Cy (z/2)  y 

kyC^  z  +  SjLs/2  +  e,  (Lj  -  2z)/2 

Solving  for  y,  and  noting  that  at  z  =  L,/2,  y  =  yy,  obtain  the  following 
equation  for  yy, 

€y  Lf 

Substituting  the  value  of  kydg  from  Eq.  22,  with  z  =  Ls/2,  Into  Eq.  26, 
and  rearranging,  the  following  quadratic  equation  Is  obtained, 

-  zjh,  +  ^  -  ri)jyy  +  Y  [sud  +  8,)  +  s/j  =  0 

For  Ny,  a:  N^i  and  Ny,  S  N^g  (tension  failure  at  both  ends)  this  can  be 
rewritten  as 
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(27) 


y?  -  2[h,  +  ^  ('3  -  *‘i)]yu  +  ^  [eu(i  +  8.)  +  s/j  =  0 

Solvliifr  for  jTy ,  the  central  deflection  required  to  crush  the  concrete 
along  the  hinge  lines  for  the  z-dlrectlon  strips  Is 


where  r^  and  rg  are  given  by  Eqs.  23a  and  24a, 

Equation  28  shows  that  Xy  depends  on  (1)  the  properties  of  the  cross 
section  at  the  ends  of  the  strip;  (2)  the  crushing  strain  of  the  concrete; 
(3)  the  span  and  thickness  of  the  slab;  (4)  the  amount  of  longitudinal 
movement  at  the  supports;  (5)  the  axial  stiffness  of  the  slab;  (6)  the 
depth  to  the  tenslcm  steel  relative  to  the  slab  thickness;  and  (7)  the 
ratio  of  f,  to  fy. 

For  certain  combinations  of  these  slab  properties,  the  quantity 
under  the  radical  In  Eq.  28  may  be  negative.  This  Indicates  Insufficient 
longitudinal  restraint  In  the  z-dlrectlon  to  develop  crushing  strain  along 
the  hinge  lines,  thus  resulting  In  an  Imaginary  value  for  y^.  In  this 
case  the  corresponding  ultimate  deflection  for  the  x-dlrectlon  (long 
direction)  strip  Is  calculated. 

Following  the  same  procedure  as  for  the  z-dlrectlon  strips,  the 
central  deflection  required  to  crush  the  concrete  along  the  hinge  lines 
for  the  x-dlrectlon  strips  Is 
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Yu 


(r^ 


■4 


N  “  **3  ^ 


ij®  -  pl^[^€u(2p  +  S^  -  2pe,  +  e,)  +  s/j 

(29) 


where  and  are  given  by  Eqs.  16a  and  17a. 

If  the  value  of  yu  from  Eq.  29  is  also  Imaginary,  the  longitudinal 
restraint  In  both  directions  is  insufficient  to  develop  the  crushing 
strain  along  the  yield  lines.  For  such  cases  the  compressive  membrane 
forces  are  insignificant,  and  the  slab  resistance  is  calculated  using 
the  yield-line  theory  previously  outlined  in  Ref.  7. 

Equations  28  and  29  do  not  apply  if  the  ultimate  thrust,  ,  is 
greater  than  the  balanced  thrust,  N^, .  However,  since  an  exact  solution 
for  such  cases  would  require  a  complicated  iterative  procedure  to  solve 
for  fg  and  along  the  yield  lines,  Eqs.  28  and  29  were  also  used  for 
those  slabs  where  Nu  >N^.  For  such  cases,  the  slabs  have  limited  ro¬ 
tational  capacity  along  the  hinge  lines,  thus  introducing  the  possibility 
of  a  brittle  collapse  mechanism  that  can  precipitate  a  premature  shear 
failure.  As  noted  previously,  this  situation  was  not  considered  in  this 
study  due  to  the  lack  of  information  regso'dlng  the  increase  in  shear 
strength  due  to  the  presence  of  compressive  membrane  forces. 

Examining  Eqs.  28  and  29  in  more  detail,  it  can  be  noted  that  y^ 
reaches  a  limiting  value  when  the  quantity  under  the  radical  equals  zero. 
In  other  words,  the  compressive  membrane  forces  induced  in  the  slab  reach 
a  maximum  value  regardless  of  the  span/thickness  ratio.  According  to 
Eqs.  28  and  29  the  upper  bound  for  yy  is  h,  +  (r3  -  r^  )/k]^  for  the  z- 
dlrectlon  strips  and  h,  +  (r4  -  for  the  x-dlrectlon  strips.  These 

upper  bounds  correspond  to  large  values  of  the  span/thickness  ratio. 
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Holwver,  examination  of  test  values  from  Refs.  9,  10,  11,  and  12  indicate 
a  limiting  value  of  much  less  than  the  slab  thickness.  Reference  10 
suggests  an  empirical  value  for  the  limiting  ultimate  deflection  of 

<  0,42  h,  .  (30) 

This  u.)pc.i  bound  was  also  used  in  this  study. 

Ultimate  Flexural  Resistance 

The  ultimate  flexural  resistance  can  now  be  determined  by  means  of 
the  work-energy  method.  The  work  done  can  be  found  by  multiplying  the 
resisting  moment  per  unit  length  normal  to  the  axis  of  rotation  by  the 
angle  of  rotation  for  each  of  the  x-directlon  and  z-directlon  strips, 
and  summing  for  all  the  strips. 

The  sum  of  the  moments  on  the  x-directlon  strips  are 
Mua  (dz)  +  Mu4  (dz)  -  Ntt,y  (dz) 
while  the  angle  of  rotation  is 


Vu 


The  work  done  for  the  x-dlrectlon  strips  is  thus  given  by 

Wx  -  +  m„4  -  N„,yj|^jdz 

where,  from  Eq.  8, 

(it)  » * « *  ek  • 


(31) 


(8) 
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The  corresponding  work  done  for  the  z-dlrectlon  strips  Is 


W 


t 


r . 

1/  yu  \ 

-  Ny.y 

l(k/2) 

dx 


where,  from  Eq.  19, 


y  = 


0 


Ls 


(32) 


(19) 


The  energy  Input  for  a  uniform  load,  q,  can  be  determined  by  multi¬ 
plying  the  total  load  acting  on  the  strip  by  the  deflection  at  the  center 
of  gravity  of  the  strip.  Thus,  the  energy  Input  for  the  x-dlrectlon 
strips  Is 


E, 


[q^  (dz)x] 


^  Qu  yu 
2  “  20L^ 


x®dz 


(33) 


where 

Similarly,  the  energy  Input  for  the  z-dlrectlon  strips  Is 


E* 


[q^  (dx)z] 


^  yu 
2  =  ~ 


z®dx 


where 


0  <  X  < 


3k 


k 

<  X  S  ^ 
2 


(35) 


(36a) 


(36b) 
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The  total  work  done  by  the  x-direction  and  z-direction  strips  is 
Obtained  by  integrating  over  the  total  area  of  the  slab.  Thus, 


Wt 


W. 


which,  by  substituting  Eqs.  31  and  32,  becomes 


Wt 


[Mug  +  Mu4 


Nuxy] 


Yu  V 
Pk/dz 


2 


[  +  My 


3 


(37) 


Substituting  the  expressions  for  Nu,  and  (Eqs.  18  and  25),  y  (Eqs. 

8  and  19),  and  x  and  z  (Eqs.  34,  36a,  and  36b)  into  Eq.  37,  and  integrating 
yields  the  following  equation  for  Wj, 
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Since  My4 ,  ,  and  Myg  are  nonlinear,  the  first  four  Integrals  are 

awkward  to  evaluate.  A  simple  expression  that  may  be  used  to  approximate 
these  Integrals  Is 


Ls/2 


(39a) 


(39b) 
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My^dx  =  (Muj  ) 


(39c) 


(M„3  )  Pk 


(39d) 


where 

Myg  =  the  ultimate  moment  resistance  about  mldthlckness 
of  sections  along  the  diagonal  produced  by  the 
average  thrust  acting  across  the  span,  [Nu, (0) 

+  N„^Ok)]/2. 

the  ultimate  moment  resistance  about  mldthlckness 
of  sections  along  the  edge  produced  by  the  average 
thrust  acting  across  the  span,  (0)  +  Nyj(pI^)]/2. 

=  the  ultimate  moment  resistance  about  mldthlckness 
of  sections  along  the  diagonal  produced  by  the 
average  thrust  acting  across  the  span,  [N^jCO) 

+  N,,(L3/2)]/2. 

Myg  =  the  ultimate  moment  resistance  about  mldthlckness 
of  sections  along  the  edge  produced  by  the  average 
thrust  actii^  across  the  span,  [^^,(0)  -f  N^,  (L,/2)]/2. 


The  fifth  and  sixth  Integrals  can  be  easily  evaluated,  since  Myj  and  My4 
are  constant  across  the  center  yield  line.  Thus, 


p)k“ui(ek) 


(39e) 


(39f) 
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where 


(01^ )  =  the  ultimate  moment  resistance  along  the 
center  produced  by  the  thrust  Nyx(Lj/2) 

-  the  ultimate  moment  resistance  along  the 
edge  produced  by  the  thrust  Ny,(L3/2). 


The  last  three  terms  of  Eq.  38  can  also  be  expressed  in  terms  of  the 
thrust  by  noting  from  Eqs,  17  and  24  that 


|[l.„.(0)  +  2I.„  (^)] 
*  2N„(Bk)] 


Sjlif 

1 

4yu 

)  ^ 

K  r 

/h 

2  L 

\  *  3 

4yu  / 

K  r 

(h 

2  L 

\  *  3 

4y«  / 

(39g) 

(rj  +  r3)j 

(39h) 

(rj  +  )j 

(391) 


Substituting  Eqs.  39a  through  391  In  Kq.  38  yields  the  following  expres¬ 
sion  for  , 


Wt  = 


j"..  *  «U4  -  [»..(»)  t  2».,<9k)]j 


4y„(l  -  23)1^ 

+  - r -  Mui(Pk)  +  »^3  0k) 


(40) 


The  total  energy  Input  can  be  determined  by  Integrating  the  energy 
Input  for  the  individual  x-dlrectlon  and  z-dlrectlon  strips  over  the 
entire  area  of  the  slab.  Thus, 
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0  0 

which,  by  substituting  Eqs.  33  and  35,  becomes 


Substituting  the  expressions  for  x  and  z  from  Eqs.  34,  36a,  and  36b,  and 
integrating,  the  total  energy  input  is  expressed  by 


Et  = 


quyu^sk 

6 


(3-20) 


(41) 


Equating  the  total  work  done  (Eq.  40)  to  the  total  energy  input  (£q.  41), 
and  solving  for  the  slab  resistance,  yields 

=  "  2n,,(0i,)]| 

+  ^  j^N„,(0)  +  2N„,  (~)]| 

+  2(1  -  2e)jM,i(eLL)  +  M,,(0Ll)  -  y„N„,  (~)jj  .  (42) 


The  coordinates  of  the  point  defining  the  ultimate  flexural  resistance 
are  now  defined  by  Eq.  28,  29,  or  30  for  y„  and  by  Eq.  42  for  q^  .  Reference 
10  found  that  the  resistance  for  values  less  than  q^  can  be  approximated 
by  the  expression 


y*  s  yu 


(43) 
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This  equation,  which  assumes  a  vertical  tangent  at  (0,0)  and  a  horizontal 
tangent  at  (q„ ,  y^),  was  determined  by  a  best  fit  of  the  experimental 
results  from  static  tests  of  square  slabs  longitudinally  restrained  on 
all  edges.  In  this  study,  this  expression  is  also  assumed  to  apply  to 
rectangular  slabs  and  to  slabs  with  partially  restrained  edges. 

Secondary  Resistance 

After  the  resistance  reaches  its  ultimate  value,  there  is  a  decrease 
in  the  resistance  as  the  deflection  increases.  Tests  show  that  the 
resistance  reaches  a  minimum  value  approximately  equal  to  the  value  ob¬ 
tained  using  yield-line  theory  (not  considering  membrane  forces).  This 
equivalence  is  reasonable  to  expect,  since  the  minimum  resistance  should 
correspond  to  a  change  of  membrane  forces  from  compression  to  tension  in 
the  central  region  of  the  slab.  The  secondary  resistance  can  thus  be 
determined  from  £q.  42  by  setting  and  equal  to  zero,  resulting 
in  the  following  equation  for  q^, 

where  and  are  the  resisting  moments  for  a  zero  thrust 

and  3  determined  from  Eq.  1. 

Expressions  previously  given  in  Ref.  7  for  the  tensile  membrane  be¬ 
havior  of  slabs  are  used  to  determine  the  deflection  at  the  intersection 
of  the  tensile  membrane  resistance  and  secondary  resistance.  These  ex¬ 
pressions,  which  were  based  primarily  on  the  development  presented  in 
Ref.  13,  assume  the  reinforcing  steel  to  act  as  a  plastic  membrane. 
Additional  assumptions,  all  of  which  are  conservative,  include: 
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(1)  The  concrete  has  cracked  throughout  its  depth,  hence 
is  incapable  of  carrying  any  load. 

(2)  All  the  reinforcing  steel  has  yielded. 

(3)  No  strain-hardening  of  the  reinforcing  steel  occurs. 

The  resulting  expression  for  the  tensile  membrane  resistance  is 


q 


y* 


(45) 


where 


g  is  the  area  (per  unit  width)  of  the  reinforcing  steel  (both  top  and 
bottom  bars)  continuous  along  the  short  span  of  the  slab.  Similarly, 

A^^  is  the  area  (per  unit  width)  of  the  reinforcing  steel  continuous 
along  the  long  span  of  the  slab.  Values  of  the  factor  are  plotted  in 
Figure  10  for  A^g/A^i^  a  1.  The  1.5  factor  was  included  to  reflect  the 
approximately  50  percent  Increase  in  test  values  of  the  tensile  membrane 
strength  over  the  calculated  values. 

Setting  q  =  q,  in  Eq.  45  and  solving  for  the  center  deflection  at 
the  intersection  of  the  secondary  resistance  and  tensile  membrane  re¬ 
sistance. 


Vt  = 


% 


k'A, 


S^dy 


(47) 


In  general,  the  secondary  deflection  is  taken  to  be  the  same  as  the 
Intersection  of  the  secondary  resistance  and  the  tensile  membrane  resistance. 


or 
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k/Ls 


SOURCE:  Rtftrenca  13. 


FIGURE  10  COEFFICIENTS  FOR  TENSILE  MEMBRANE  RESISTANCE 
IN  TWO-WAY  ACTION  SLABS 
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y.  =  vt 


(48) 


q,L| 


S  y 


Examination  of  test  results  Indicate  that  for  certain  slabs,  primarily 
those  with  low  tensile  membrane  resistance,  the  value  given  by  Eq.  48  Is 
much  larger  than  the  actual  secondary  deflection.  Based  on  a  study  of 
test  results,  an  upper  bound  of 


y,  s  3  y„ 


(49) 


Is  placed  on  the  secondary  deflection. 

The  coordinates  of  the  point  defining  the  secondary  resistance  are 
now  defined  by  Eq.  44  for  q,  and  by  Eq.  48  or  49  for  y, .  In  Ref.  10 
resistance  for  values  between  y^  and  y,  was  approximated  by  the  expression 


q  = 


5.  [i  r.  .  (Ml  j 

2  I  Uu/  L  v‘>u/J  L  y.  “  yu  J) 


(50) 


This  expression  was  also  used  In  this  study  to  describe  the  resistance- 
deflection  behavior  for  this  region. 

For  the  special  case  of  a  slab  with  no  reinforcement,  the  secondary 
resistance  Is  zero.  For  this  case  the  secondary  deflection  Is  defined 
to  be  equal  to  the  slab  thickness. 


y,  =  h,  .  (51) 

The  resistance  for  values  between  qy  and  q.  Is  still  given  by  Eq.  50. 


Tensile  Membrane  Resistance  and  Collapse 

After  the  deflection  y^  Is  reached,  the  tensile  membrane  region 
begins  to  grow  outward  towards  the  supports.  Cracks  begin  to  penetrate 
the  entire  depth  of  the  concrete,  and  yielding  of  the  reinforcement 
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spreads  throughout  the  slab.  The  top  and  bottom  reinforcement  acts  as  a 
tensile  membrane,  causing  the  slab  resistance  to  Increase  with  the  center 
deflection,  as  shown  by  Eq.  42.  This  continues  until  the  reinforcement 
ruptures.  Reference  13  reported  that  a  safe  maximum  value  for  the  central 
deflection  Is  0.10  of  the  short  span.  Examination  of  test  results  Indi¬ 
cates  that  a  more  realistic  value  for  collapse  of  the  slab  Is  approximately 
0.15  of  the  short  span.  Thus,  for  the  present  study,  the  failure  deflec¬ 
tion  Is  taken,  as 

y,^  =  0.15  Lj  .  (52) 

For  slabs  unrestrained  In  the  short  direction,  Eq.  52  does  not  apply. 

For  these  cases,  the  failure  deflection  Is  taken  as  0.15  of  the  long  span, 

=  0.15  •  (53) 


Summary 

The  resistance  function  for  reinforced  concrete  slabs  with  longi¬ 
tudinal  edge  restraint,  as  previously  Illustrated  In  Figure  1,  with  the 
corresponding  equations  is  summarized  as  follows. 

Compressive  Membrane  Resistance  (0  S  y^  ^  y„  ) 


bu  —  Eq-  42 

y^  -  Eq.  28,  Eq.  29,  or  Eq.  30 
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Secondary  Resistance  (y^  i  ^  yt^ 


(50) 


q,  -  Eq.  44 

y^  -  Eq.  48,  Eq.  49,  or  Eq.  51 


Tensile  Membrane  Resistance  (y^  S  y,  ^  y^^) 


q  =  K 


yo 


(45) 


k/  -  Eq.  46 

y,^  Eq,  52  or  Eq.  53 

After  develoimient  of  the  resistance  function  for  longitudinally  re¬ 
strained  reinforced  crmcrete  slabs,  the  computer  programs  prepared  under 
the  previous  study  (Ref.  7)  were  modified  accordingly.  This  required 
transforming  the  slab  Into  an  equivalent  slngle-degree-of-freedom  dynamic 
system  by  the  use  of  transformation  factors  for  the  load,  resistance, 
and  mass.  The  equation  of  motion  is  solved  cm  a  computer  using  numerical 
integration  procedures.  Analytical  predictions  of  the  static  and  dynamic 
response  of  restrained  slabs  were  then  generated  for  comparison  with 
various  published  experimental  data.  The  results  of  the  analytical  and 
experimental  correlations  are  presented  In  the  Appendix.  In  general,  the 
correlations  Justify  the  use  of  the  compressive  membrane  mode  for  pre¬ 
dicting  the  collapse  of  floor  systems  In  existing  structures  that  are 
longitudinally  restrained  at  the  edges,  such  as  the  interior  panels  of  a 
multibay  floor  system. 
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Ill  WOOD-JOIST  FLOOR  SYSTEMS 


Introduction 

Wood-Joist  floor  systems  are  complex  stnictural  systems  that  are 
not  readily  amenable  to  precise  theoretical  analysis  because  of  the  wide 
variability  in  the  mechanical  properties  of  wood,  even  within  the  same 
species  and  grade,  and  the  variability  in  construction  methods  and  details. 
Even  though  it  is  known  that  the  finish  flooring  and  subflooring  con¬ 
tribute  to  the  strength  and  stiffness  of  floor  systems,  design  codes 
invariably  specify  maximum  stresses  and  deflections  for  the  design  of 
the  Joists  to  resist  the  total  floor  load  without  regard  to  other  factors. 
The  allowable  stress  specified  in  the  design  codes  are  quite  low,  since 
95  percent  of  the  individual  load-carrying  members  in  a  given  lumber  grade 
can  safely  support  the  design  load  (Ref.  14).  For  example.  Ref.  15 
specifies  an  allowable  unit  stress  in  the  extreme  fiber  in  bending  of 
2050  psi  for  the  highest  grade  Douglas  fir,  coast  region  Joist  and  plank 
stress-graded  lumber,  whereas  Ref.  16  indicates  a  modulus  of  rupture  in 
excess  of  12,000  psi  for  small  wood  beam  specimens  in  static  bending 
tests  of  the  same  species. 

As  noted  in  various  research  publications,  design  codes  do  not  take 
into  account  the  effect  of  the  Interaction  between  individual  members  of 
a  completed  floor  system,  and  the  specifying  of  relatively  low  allowable 
stresses  for  Joist  design  is  equivalent  to  requiring  floor  system  design 
to  be  based  on  the  weakest  link.  On  the  other  hand,  the  available  test 
Information  shows  that  the  strength  of  a  floor  system,  beyond  that  of  the 
weakest  Joist,  is  largely  dependent  on  the  interaction  effects  between 
floor  elements.  For  example,  the  finish  and  subflooring  both  contribute 
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to  the  strength  and  stiffness  of  the  Joist  through  composite  or  T-beam 
action.  However,  the  primary  effect  of  the  flooring  is  a  result  of  its 
contribution  to  load  sharing  between  individual  Joists.  The  wood  in¬ 
dustry's  concept  of  load  sharing  has  been  described  in  Ref.  14  as  being 
a  reduction  in  the  effective  variability  between  individual  joists  when 
they  are  grouped  into  a  structural  system,  an  increase  in  the  effective 
strength  provided  by  the  mutual  constraint  resulting  from  Joining  members 
into  an  indeterminate  structure,  and  an  Increase  in  the  effective  strength 
by  the  local  reinforcement  of  defects  by  adjacent  members. 

Although  design  procedures  do  not  permit  assigning  a  value  to  load 
sharing,  as  mentioned,  the  enhancement  of  the  floor  system  strength  as 
a  result  of  the  interacting  effects  among  the  various  elements  in  a  wood- 
joist  floor  system  has  been  conclusively  demonstrated  in  many  tests  of 
actual  floors  (e.g..  Refs.  17  through  22).  The  intent  of  this  discussion 
is  not  to  Judge  the  merits  of  any  code  requirements,  but  merely  to  point 
out  that  for  predicting  the  collapse  of  wood-joist  floors,  it  is  necessary 
to  include  effects  that  are  not  treated  directly  in  any  design  procedure. 

The  development  of  a  realistic  resistance  function  to  represent  the 
response  of  a  wood-Jolst  floor  system  is  complicated  by  the  large  vari¬ 
ability  of  material  properties  within  any  pcorticular  wood  species  (even 
among  similar  wood  specimens  specifically  selected  for  tests),  and  by 
the  highly  complex  indeterminate  structural  action  that  occurs  between 
various  elements  of  the  floor  system.  Although  the  limited  theoretical 
studies  of  the  composite  floor  behavior,  such  as  that  presented  in  Ref. 

23  to  model  the  slip  in  nailed  plate-rib  Joints,  have  been  correlated 
with  specific  test  data,  the  unknowns  involved  in  the  evaluation  of  wood- 
joist  floors  in  existing  buildings  do  not  warrant  the  development  of  a 
sophisticated  mathematical  model.  Instead,  the  approach  in  this  study 
was  to  examine  sufficient  analytical  and  experimental  information  to 
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permit  the  development  of  a  relatively  simplified  model  that  would  ade¬ 
quately  reflect  the  complex  composite  action  of  the  floor  system. 

A  wood  beam,  or  Joist,  composed  of  perfectly  clear  wood  without 
defects,  and  subjected  to  pure  bending  forces  as  a  result  of  an  Increasing 
uniform  static  load,  goes  through  two  distinct  phases  prior  to  collapse. 

The  first  phase  Is  elastic,  where  a  plane  section  remains  plane,  and  a 
linear  relationship  exists  between  stress  and  strain  throughout  the  beam 
cross  section.  As  the  stress  level  Is  Increased,  the  extreme  fibers  on 
the  compression  side  reach  their  proportional  limit.  The  second  phase  Is 
therefore  characterized  by  Inelastic  behavior  on  the  compression  side,  and 
a  linear  elastic  behavior  on  the  tension  side.  Typical  stress-strain 
curves  for  wood  in  tension  and  compression  are  shown  on  Figure  11;  note 
that  the  tension  portloi?  of  the  stress-strain  relationship  remains  essen¬ 
tially  linear  up  to  failure.  Since  the  ultimate  tensile  strength  of 
clear  wood  Is  much  greater  than  its  ultimate  compressive  strength,  a 
redistribution  of  stress  occurs  across  the  beam  section  as  the  bending 
moment  increases  beyond  the  proportional  limit  In  compression.  To  satisfy 
the  conditions  of  eq\iillbrlum  of  internal  forces  at  a  section,  the  neutral 
axis  shifts  towards  the  tension  side  (Refs.  24  and  25).  Test  data  also 
Indicate  that  a  plane  section  of  a  wood  beam  In  pure  bending  remains 
approximately  plane  up  to  tensile  failure. 

As  noted  In  Ref.  25,  the  second  degree  parabola  has  been  found  to  be 
the  best  approximation  to  the  compressive  stress-strain  distribution  above 
the  proportional  limit.  The  theoretical  stress  distribution  In  a  wood 
beam  near  failure  would  therefore  appear  as  shown  In  Figure  12. 

Resistance  Function 

As  discussed  above,  a  wood-jolst  floor  system  is  a  complex  structural 
system  because  of  both  the  effects  of  Interaction  between  the  floor  elements, 
and  the  wide  variability  In  the  properties  of  wood.  In  addition. 
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FIGURE  11  STRESS-STRAIN  CURVES  FOR  AXIAL  TENSION 
AND  COMPRESSION  IN  WOOD 


COMPRESSION 


FIGURE  12  STRESS  DISTRIBUTION  NEAR  FAILURE  FOR 
RECTANGULAR  WOOD  BEAMS  IN  FLEXURE 
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predicting  the  collapse  of  existing  floor  systems  under  blast  loading  is 
further  complicated  by  the  variability  in  construction  techniques, 
especially  in  the  relatively  unsupervised  home  construction  Industry, 
and  by  the  fact  that  the  properties  of  the  wood  used  in  any  analysis  are 
usually  published  average  values  for  a  wood  species  rather  than  actual 
values  based  on  knowledge  of  the  specific  wood  used.  The  differences 
between  the  published  pro'^ertles  of  wood  and  those  in  a  specific  existing 
floor  system  would  be  expected  to  show  a  much  greater  variation  than  for 
structural  steel,  or  even  for  concrete.  However,  the  use  of  published 
average  property  values  for  wood  may  be  adequate  because  of  the  averaging 
effect  resulting  from  the  interaction  among  the  elements  of  a  floor 
system. 

In  any  event,  as  stated,  the  unknowns  involved  in  calculating  the 
dynamic  response  and  collapse  of  wood-joist  floor  systems  do  not  Justify 
the  development  of  a  complex  mathematical  model.  The  approach  therefore 
was  to  use  established  theoretical  procedures,  and,  where  feasible,  to 
modify  the  procedures  to  reflect  the  results  of  tests  of  actual  floor 
systems.  The  method  adopted  was  similar  to  that  used  for  wall  elements 
(Refs.  1  and  4),  and  for  reinforced  concrete  floor  systems  (Ref.  7). 

That  is,  the  resistance  function  was  established  for  wood-joist  floor 
systems  by  assuming  that  the  deflected  shape  of  the  member  under  dynamic 
load  is  identical  to  that  under  a  uniform  static  lo«d  and  that  the  distri¬ 
bution  of  the  restoring  force  and  dynamic  load  is  whe  same.  The  member 
is  then  transformed  into  an  equivalent  slngle-degree-of-freedom  dynamic 
system  by  the  use  of  transformation  factors  for  the  load,  resistance,  and 
mass.  The  equation  of  motion  is  solved  on  a  ccxnputer  using  a  numerical 
integration  procedure. 

The  resistance  function  for  wood-joist  floor  systems  was  assumed  to 
consist  of  elastic  and  inelastic  phases  as  shown  on  Figure  13,  and  was 
determined  only  for  the  case  of  a  simply  supported,  uniformly  loaded  floor. 
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RESISTANCE 


FIGURE  13  RESISTANCE  FUNCTION  FOR  A  SIMPLY  SUPPORTED 
WOOD-JOIST  FLOOR  SYSTEM 
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Although  scabs  are  sometimes  used  to  reinforce  wood  Joists  at  Interior 
supports  (e.g.,  Ref.  26),  Joists  are  not  generally  continuous  over  the 
support,  and  the  scab  cannot  be  expected  to  develop  the  bending  strength 
of  the  Joist.  Also,  although  tests  Indicate  that  a  celling  on  the  bottom 
of  a  wood-Jolst  floor  has  a  stiffening  effect  (Refs.  26  and  27),  a 
celling  was  not  Included  in  the  mathematical  model  since  the  primary 
Interest  was  in  floors  over  basement  areas  where  ceilings  are  usually 
not  used. 

Elastic  Phase 

The  maximum  elastic  resistance,  ,  of  a  simply  supported  wood-Jolst 
floor  with  a  uniformly  distributed  load  Is  developed  when  the  moment  at 
the  center  section  Is  a  maximum,  or 

q^sL® 

M,  =  ,  (54) 

where  s  Is  the  Joist  spacing. 

Since  a  linear  relationship  Is  assumed  to  exist  between  the  stress 
and  strain  during  the  Initial  elastic  phase,  the  extreme  fiber  stress  Is 
equal  to 

Mhj  /2 

f,  = 

Test  data  (e.g..  Ref.  28)  have  Indicated  that  the  effective  moment  of 
Inertia,  I,,  of  a  wood-joist  floor  system  Is  greater  than  the  moment  of 
Inertia  of  the  Joists,  ,  but  less  than  that  obtained  from  a  T-section 
(Figure  14)  if  full  composite  action  between  the  Joist  and  flooring  Is 
assumed.  For  this  study,  the  effect  of  the  composite  action  Is  accounted 
for  in  the  elastic  phase  by  assuming  that  the  effective  moment  of  Inertia 
Is  equal  to  the  moment  of  Inertia  of  the  Joist  times  a  coefficient  that 
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FIGURE  14  CROSS  SECTION  OF 


/ 


FLOORING 
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Is  determined  from  the  published  test  data  on  the  response  of  wood-Jolst 
floors  In  the  elastic  range,  or 

I,  =c,Ij  ,  (56) 


or,  for  a  rectangular  Joist, 


I. 


(57) 


Substituting  Eqs.  55  and  57  Into  Eq.  54,  and  rearranging  terms,  the 
maximum  elastic  resistance  for  a  wood-Jolst  floor  system  Is 


Qi  = 


4f,c,bh^ 

3sL® 


(58) 


The  maximum  deflection  for  the  elastic  phase  Is 
5qisL‘ 

'  32EjC,bh^ 


(59) 


The  problem,  of  course.  Is  to  determine  the  point  In  the  response 
of  the  floor  system  where  It  Is  no  longer  linearly  elastic;  this  Is 
primarily  a  problem  of  determining  a  value  for  the  extreme  fiber  stress, 
f^,  for  use  In  Eq.  58  that  reflects  the  upper  limit  of  the  elastic  phase. 
Initially,  In  this  study  It  was  assumed  that  the  elastic  phase  existed 
until  the  extreme  fiber  stress  In  compression  reached  the  proportional 
limit,  that  Is  where  f^  =  f^gp  In  Eq.  58.  However,  comparisons  of  the 
analytical  predictions  with  experimental  data  (Refs.  21  and  29)  Indicated 
that  a  better  correlation  could  be  obtained  by  estimating  the  failure 
strength  of  the  weakest  Joist  (first  break),  and  using  this  value  as  the 
maximum  elastic  resistance  of  the  floor  system.  In  fact,  test  data  In 
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Ref.  21  show  this  to  be  the  case,  at  least  for  the  lower  grades  of  wood* 

♦ 

Joists. 

As  a  result  of  load  sharing,  the  deflected  shape  of  a  floor  system, 
which  consists  of  a  number  of  Joists  connected  by  flooring,  is  about  the 
same  throughout  the  floor  when  uniformly  loaded  (Refs.  28  and  30).  The 
deflected  shape  of  the  weakest  Joist  will  approximate  that  of  its  stronger 
neighbors,  even  though  it  may  have  much  less  bending  stiffness,  and 
therefore  carry  less  load.  For  any  specific  wood  species,  the  strength 
of  a  Joist  in  bending  is  influenced  by  a  number  of  factors,  such  as  size 
and  location  of  knots,  shakes  and  splits,  moisture  content,  density,  and 
slope  of  the  grain  (Refs.  16,  31,  and  32).  To  determine  the  maximum 
elastic  resistance,  it  was  assumed  that  the  bending  strength  of  the 
weakest  Joist  was  limited  by  a  defect  resulting  from  the  largest  per¬ 
missible  knot^  located  at  the  bottom  edge  of  the  wide  face  in  the  middle 
third  of  the  Joist  span  as  specified  in  the  appropriate  standard  grading 
rules  for  the  wood  grade  and  species  used  (e.g..  Refs.  32  and  33). 

The  maximum  elastic  resistance  at  first  break  is  calculated  by  using 
Eq.  58,  and  substituting  f„  ~  f,  and  Cf  -  1.0,  or 

'**.'’‘'1  .  (58a) 

'  Sat," 

The  modulus  of  rupture,  f , ,  is  that  obtained  from  standard  sources  for 
small  clear  specimens  (Ref.  34),  and  adjusted  by  the  appropriate  factors 
for  knots,  seasoning  or  moisture  content,  and  size  effect  from  Ref.  31. 

♦ 

Limited  data  on  the  collapse  of  floors  constructed  of  a  higher  grade 
Douglas  fir  in  Refs.  18  through  20  also  indicate  that  first  break  is  a 
good  measure  of  the  maximum  elastic  resistance. 

^As  noted  in  Ref.  16,  for  design  purposes  the  reduction  in  strength  for 
vcurious  defects  is  not  considered  as  cumulative.  Because  of  the  dis¬ 
tribution  of  defects,  this  is  also  a  reasonable  assumption  for  the  cal¬ 
culation  of  the  floor  resistance  at  first  break. 
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Inelastic  Phase 


To  develop  the  equation  for  the  resistance  function  during  the  in¬ 
elastic  or  plastic  phase  of  the  response  of  wood-joist  floor  systems,  it 
was  assumed  that  point  2  on  Figure  13  represented  the  maximum  theoretical 
resistance  for  a  clear  wood  specimen.  The  resistance  was  assumed  to  vary 
linearly  between  the  maximum  elastic  and  inelastic  resistances,  and  there¬ 
fore,  it  was  only  necessary  to  determine  the  value  of  the  maximum  inelastic 
resistance. 

Studies  have  shown  that  as  the  load  on  a  wood  beam  is  increased,  a 
local  failure  first  occurs  in  the  compression  zone,  and  is  accompanied  by 
a  shifting  of  the  neutral  axis  towards  the  tension  side  (e.g.,  Refs.  24, 

25,  and  35).  For  clear  specimens  without  defects,  the  collapse  of  the 
beam  results  from  a  tensile  fracture.  Although  several  theories  have  been 
postulated  to  describe  the  stress  distribution  across  a  wood  beam  at 
failure,  the  second  degree  parabola  as  presented  in  Ref.  25,  and  shown 
on  Figure  15,  was  adopted  in  this  study. 

In  the  development  of  the  equation  for  the  ultimate  plastic  moment 
at  failure  of  a  wood  beam,  the  following  assumptions  were  made  in  Ref.  25: 

•  Cross  sections  of  the  beam  remain  plane  up  to  failure. 

•  The  tensile  stress-strain  relailonship  is  linear. 

•  The  compressive  stress-strain  relationship  is  linear 

up  to  the  proportional  limit,  and  thereafter  is  approximated 
by  a  second  degree  parabola. 

•  The  downward  shift  in  the  neutral  axis  at  failure  can  be 
described  by  the  relationship  of  the  ultimate  compressive 
and  tensile  stress  distributions. 

The  ultimate  internal  resisting  moment  for  a  rectangular  wood  beam 
in  pure  bending  can  be  formulated  from  the  stress-strain  relationship  of 
wood  in  tension  and  compression  and  the  general  equations  of  static 
equilibrium  across  a  section,  where  the  sum  of  the  normal  stresses  is 
equal  to  zero,  or 
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FIGURE  15  THEORETICAL  STRESS  DISTRIBUTION  IN  CLEAR  WOOD 
BEAM  SPECIMEN  AT  TENSILE  FRACTURE 
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f,dA 


0 


and  the  moment  of  the  normal  stress  about  the  neutral  axis  Is  equal  to 
the  bending  moment  M,  or 


I 


A 


ydA 


M 


For  the  assumed  case  of  a  rectangular  beam  with  the  compression  zone 
described  by  a  second  degree  parabola  and  the  tension  zone  by  a  straight 
line  as  shown  In  Figure  15,  the  equations  of  static  equilibrium  become 


oh 


and 


-  — 

b  I  f,dy  =  0 

•'-hj  (l-o) 


f,dy  .  M 

^  •_  Mm.  _  V 


-hj (1-0) 


(60) 


(61) 


At  the  development  of  the  ultimate  bending  resistance  of  the  beam,  the 
best  fit  to  the  test  data  for  the  distribution  of  the  stress  In  the  com¬ 
pressive  zone,  from  Ref.  25,  Is  a  second  degree  parabola  of  the  form 


o-h? 


- 


2f 


wou 


oth, 


(62) 
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and  the  distribution  of  the  tensile  stress  is  a  straight  line  of  the  form 


X  =  - 


‘■(•to 


hj (1-a) 


(63) 


Since,  from  Eq.  60,  the  sum  of  the  areas  of  the  tensile  and  com¬ 
pressive  stresses  across  a  section  are  equal  to  zero,  then 


F  =  b 


0 

•1 


y  dy  =  0 


-hj (l-o) 


and  by  integrating  and  collecting  terms. 


3f 


a  = 


wt  u 


3f j  +  4f^gy 


For  n  —  , 


a 


3n 

3n  4  4 


(64) 


The  ultimate  internal  resisting  moment  of  the  stress  diagram  in 
Figure  15  by  Eq.  61  is 


(1-a) 


hj  (1-a) 


y®  dy 


f 
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and  by  integrating  and  collecting  terms 


(65) 


Since  O'  =  3n/(3n  +  4),  then 


nbh?f„au 

^  - i_2»« —  (45„  +  64) 

12(3n  +  4)3 


(66) 


To  determine  the  ultimate 
it  was  assumed  that  the  effect 
floor  system  was  negligible  at 
ultimate  plastic  resistemce  of 
spacing  s  is 


resistance  of  a  wood-joist  floor  system, 
of  the  flooring  on  the  resistance  of  the 
collapse  of  the  Joists.  Therefore,  the 
a  wood-joist  floor  system  with  Joist 


‘lu 


8Mu 

sl^ 


(67) 


By  substituting  Eq.  66  into  Eq.  67,  the  ultimate  resistance  becomes 


% 


2nbh3f,eu 
3sL3(3n  +  4)S 


(45n  +  64) 


(68) 


The  ultimate  deflection,  y^,  at  the  development  of  the  ultimate 
bending  resistance  of  a  wood-joist  floor  is  found  from  the  usual  equation 
for  the  deflection  of  a  beam 


Sq^sL* 
"  384 (El), 


(69) 


However,  since  the  quantity  (EI),^,  at  large  Inelastic  deflections  of  a 
wood  beam  is  not  equal  to  the  flexural  rigidity,  E^l^,  it  is  necessary  to 
develop  a  rational  basis  for  calculating  (El),^^. 
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For  a  beam  subjected  to  pure  bending  In  the  elastic  region,  the 


relationship  between  the  radius  of  curvature,  p,  and  the  moment  Is 


M.Si 

P 


and  In  the  Inelastic  region  can  be  described  (Ref.  36)  as 


M  = 


(El)' 


(70) 


It  can  also  be  determined  that  the  radius  of  curvature 


h,(l  -  a) 

p  =  -i— -  , 


(71) 


and  the  strain 


c  = 


'■wt 


(72) 


By  substituting  Eqs.  71  and  72  Into  Eq.  70,  and  rearranging  terms 


h,(l  -  a)EM 
(El)'  =  -i— ; - 


At  the  development  of  the  ultimate  moment  at  the  center  sectl<m  of  a 
wood  beam,  f^^  s  f^tu  =  nf^gy,  M  =  My,  and  since  E  =  Ej  then 


(El)'  = 


bh^Ej 


3(3n  +  4)3 


(45n  +  64) 


(73) 


For  a  simply  supported  wood  beam  at  ultimate  deflection,  the  quantity 
(El) '  applies  for  the  center  section  only,  while  the  elastic  Ej Ij  applies 
at  the  supports.  In  this  study,  the  quantity  (El)^^i  used  to  determine 
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the  deflection  at  development  of  the  ultimate  resistance  In  a  wood  beam 
is  assumed  to  be  the  average  £1  at  the  center  and  supports,  or 


(El) 


•  Vg 


Ejlj  +  (El)' 
2 


Since,  for  a  rectangular  section 

bh^ 

=  li  ' 


then  by  substituting  Eqs.  73  and  75  into  Eq.  74 


(El) 


a  V  K 


4(45n  +  64)' 
24  (3n  +  4)3  _ 


(74) 


(75) 


(76) 


The  use  of  £q.  69  to  determine  the  deflection  of  a  wood  beam  at  the 
development  of  the  ultimate  plastic  resistance  of  the  beam  tacitly 
assumes  that  the  shape  of  the  elastic  curve  for  simply  supported  beams 
is  valid  up  to  collapse.  This  assumption  is  necessary,  since  none  of  the 
publications  reviewed  presented  information  on  the  deflected  shape  of 
wood  beams  at  or  near  collapse.  In  any  event,  since  Eq.  69  is  used  only 
to  estimate  the  ultimate  deflection  and  not  to  predict  collapse  (as  noted 
in  the  next  subsection),  the  difference  in  the  assumed  and  actual  deflected 
shapes  is  not  felt  to  be  of  great  importance. 

Also,  the  above  equations  consider  only  the  bending  resistance  of  a 
wood-joist  floor  system,  although  a  horizontal  shear  failure  can  limit 
or  modify  the  resista;xe  function.  Studies  of  the  effect  of  the  span-to- 
depth  ratio  indicate  that  failure  in  horizontal  shecu*  is  not  likely  to 
occur  prior  to  a  tensile  failure  in  bending  for  the  span-to-depth  ratios 
usually  found  in  actual  wood-joist  floor  systems  (Ref.  37).  For  example, 
the  test  results  for  full-scale  floors  in  Refs.  17  through  20  show  that 
all  12  of  the  floors  tested  failed  in  tension,  or  rupture,  of  the  Joists 
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(frequently  Influenced  by  knots),  and  only  a  few  Joists  experienced  a 
horizontal  shear  failure  in  addition  to  tension  failure. 

Failure  Criterion 

As  mentioned  in  the  introduction!  wood-joist  floor  systems  are  con¬ 
structed  of  Joists  with  laT^e  differences  in  strength  and  stiffness 
characteristics,  even  wdien  the  Joists  are  specifically  selected  for  uni¬ 
formity  as  test  specimens.  For  example,  compression  and  tensile  strength 
values  obtained  from  tests  of  a  number  of  small,  clear  specimens  of  a 
single  wood  species  exhibit  a  wide  variation,  and  generally  conform  to 
a  typical  normal  type  of  frequency  distribution  (e.g..  Ref.  25).  However, 
regardless  of  their  Individual  stiffnesses,  when  Joists  are  part  of  a 
floor  system,  test  results  Indicate  that  the  deflection  of  the  Joists  are 
approximately  the  same  under  uniform  load  conditions  (e.g..  Refs.  21  and 
30);  the  effect  of  load  sharing  Is  an  effective  reduction  In  the  varia¬ 
bility  among  Individual  Joists.  When  tested  as  Individual  beams,  many 
commercially  available  wood  Joists  apparently  fall  as  a  result  of  a 
local  defect  while  still  within  the  elastic  range  (Ref.  14);  therefore, 
load  sharing  also  has  the  effect  of  Increasing  the  strength  of  a  floor 
system  over  that  of  the  weakest  Joist.  However,  because  of  local  defects 
and  the  variability  of  strength  among  Joists,  the  ultimate  strength  of  a 
floor  system  would  represent  some  average  value  rather  than  approach 
either  the  strength  of  the  weakest  Joist  or  the  theoretical  strength  of 
clear  specimens.  That  Is,  although  the  validity  of  Eq.  68  for  predicting 
the  ultimate  static  resistance  has  been  experimentally  verified  for  clear 
wood  specimens,  the  ultimate  resistance  of  clear  specimens  Is  not  expected 
to  represent  the  strength  of  wood-joist  floors  in  actual  buildings, 
especially  for  home  construction  where  utility  or  standard  grade  lumber 
Is  generally  used. 
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The  development  of  a  rational  failure  criterion  for  wood-Jolst 
floors  Is  complicated  by  a  number  of  factors,  some  of  which  are  mentioned 
above.  Also,  because  of  the  large  number  of  variables  Involved  In  the 
response  and  collapse  of  a  floor  system,  the  available  test  data  are 
Insufficient  to  establish  a  statistical  basis  for  developing  a  failure 
criterion.  Instead,  for  this  study  It  was  necessary  to  rather  arbitrarily 
establish  a  failure  criterion  without  the  benefit  of  adequate  analytical 
or  experimental  Information. 

Therefore,  to  predict  the  collapse  of  wcod-Jolst  floor  systems  In 
the  program  for  evaluation  of  existing  structures.  It  was  assumed  that 
the  resistance  at  collapse  would  occur  between  the  maximum  elastic  and 
Inelastic  resistances  (Points  1  and  2  In  Figure  13,  respectively).  Since 
no  rational  method  was  available  for  calculating  deterministically  the 
maximum  resistance  developed  In  a  floor  system  at  collapse,  the  collapse 
resistance  between  Points  1  and  2  was  described  In  the  form  of  a 
probability  distribution.  An  examination  of  the  experimental  data  In¬ 
dicated  that  the  maximum  resistance  developed  by  a  floor  system  at 
failure  would  be  closer  to  the  maximum  elastic  resistance  (Eq.  58a)  than 
It  would  be  to  the  ultimate  plastic  resistance  (Eq.  68);  this  suggests 
that  the  probability  distribution  for  the  failure  resistance  would  tend 
to  be  skewed  towards  the  maximum  elastic  resistance.  Although  the  0- 
dlstrlbutlon  Is  well  suited  for  describing  a  skewed  distribution  with 
finite  end  points  (Ref.  38),  only  the  normally  distributed  probability 
furctlon  was  readily  available  for  use  in  this  study.  Therefore,  it  was 
assumed  that  the  probability  of  occurrence  of  the  failure  resistance  for 
wood-Jolst  floor  systems  was  normally  distributed  between  Points  1  and  2 
on  Figure  13,  and  the  skewness  was  accounted  for  Indirectly  as  follows. 

The  possibility  is  negligible  that  the  maximum  resistance  at  failure 
for  an  actual  floor  system  would  be  greater  than  or  equal  to  the  theoreti¬ 
cal  ultimate  plastic  resistance  for  clear  wood  specimens  as  expressed 
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by  Eq.  68;  therefore,  the  probability  that  the  maximum  resistance  would 
equal  the  ultimate  resistance  was  set  very  high,  say  0.999.  On  the 
other  hand,  the  possibility  that  the  maximum  resistance  of  a  floor  would 
be  less  than  the  elastic  resistance  at  the  theoretical  first  break  of  the 
weakest  Joist,  as  expressed  by  Eq.  58a,  would  be  relatively  low;  there¬ 
fore,  for  this  case,  the  probability  of  occurrence  was  set  at  0.05  or 
0.10.  Furthermore,  in  the  dynamic  analysis  of  a  specific  floor  system, 
values  of  the  floor  resistance  above  or  below  the  limits  selected  were 
discarded.  The  truncated  normal  density  function  used  for  estimating 
the  maximum  resistance  at  failure  of  wood-joist  floor  systems  is  shown 
on  Figure  16. 

To  determine  the  probability  distribution  for  the  incipient  collapse 
overpressure  for  wood-Jolst  floor  systems,  Monte  Carlo,  or  simulation 
techniques  as  described  in  Ref.  4  for  wall  elements  was  used.  Briefly, 
the  technique  uses  a  set  of  mathematically  simulated  floors,  each  of 
which  possesses  the  characteristics  of  some  real  floor  to  determine  an 
approximate  distribution  of  the  incipient  collapse  overpressure.  The 
set  of  simulated  floors  is  prepared  by  selecting  the  parameters  to  be 
varied,  such  as  the  maximum  resistance  at  collapse,  and  determining  the 
\  alues  of  these  parameters  by  randomly  sampling  their  corresponding 
probability  distribution  functions.  Each  simulated  floor  is  then  analyzed 
dynamically  using  the  randomly  selected  values  and  the  deterministic 
equations  developed  previously.  The  results  of  the  analyses  of  the  set 
of  floors  provide  a  probability  distribution  of  the  Incipient  collapse 
overpressure. 
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IV  DYNAMIC  ANALYSIS  OF  NORTH  CAROLINA 
NATIONAL  BANK  BUILDING 


Introduction 

As  noted  In  Section  I,  a  continuing  concern  In  evaluating  the  col¬ 
lapse  overpressure  of  existing  buildings  has  been  the  relative  blast 
strength  of  the  exterior  walls  and  frames  of  multistory  buildings.  To 
predict  the  collapse  overpressure  of  the  exterior  walls  for  the  existing 
NFSS  buildings  analyzed  in  this  program  (Refs.  5  and  6),  it  was  assumed 
that  the  structural  frame  did  not  collapse  at  a  lower  overpressure  than 
that  predicted  for  the  exterior  wall.  For  weak-walled  buildings,  such 
an  assumption  is  reasonable.  In  fact,  It  Is  often  assumed  for  the  anal¬ 
ysis  of  blast-loaded  frame  buildings  that  the  exterior  walls  can  be  con¬ 
sidered  as  frangible,  and  therefore,  that  the  wall  loading  transferred 
to  the  frame  can  be  approximated  by  an  Impulse  loading.  However,  for 
many  of  the  actual  buildings  analyzed,  the  strength  of  the  exterior 
walls  under  blast  loading  was  sufficiently  high  to  make  It  doubtful  that 
the  frame  could  survive  at  the  overpressure  level  required  to  collapse 
the  walls.  For  example,  at  the  present  time  SRI  has  analyzed  a  total  of 
137  exterior  wall  cases  for  59  NFSS  buildings.  As  shown  on  Figure  17, 
the  predicted  Incipient  collapse  overpressure  for  these  walls  ranged 
from  less  than  1  psi  to  over  40  psi,  with  50  percent  of  the  walls  pre¬ 
dicted  to  collapse  at  an  incident  overpressure  level  greater  than  6  psi. 
The  strength  of  the  exterior  walls  is  Important  in  calculating  the  col¬ 
lapse  of  the  frame,  since  for  a  given  overpressure  level,  the  blast 
loading  on  the  total  area  of  a  nonfailing  wall  can  be  much  more  severe 
than  the  blast  loading  on  the  frame  alone  plus  an  impulse  loading  from 
a  frangible-type  wall. 
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FIGURE  17  CUMULATIVE  FREQUENCY  DISTRIBUTION  OF  THE  MEAN  VALUES 
OF  THE  COLLAPSE  OVERPRESSURES  FOR  THE  EXTERIOR  WALLS 
OF  59  BUILDINGS 
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To  Investigate  the  relative  strength  of  the  exterior  walls  and  frame 
of  a  building  would  require  a  comprehensive  computer  program  that  Included 
Inelastic  response  of  the  frame  under  dynamic  loading  and  realistic  frame 


collapse  mechanisms.  Since  such  a  program  was  not  available  for  this 
study,  an  existing  computer  program  for  analyzing  the  elastic  and  Inelas¬ 
tic  dynamic  response  of  two-dimensional  structural  frames  (Ref.  8)  was 
used  to  analyze  the  frame  of  the  North  Carolina  National  Bank  Building, 
This  building  was  selected  for  study  since  a  previous  analysis  of  the 
exterior  walls  (Ref.  6)  had  Indicated  that  It  was  a  strong-walled  building. 

Building  Description 

The  North  Carolina  National  Bank,  constructed  In  1922,  Is  located 
on  South  Main  Street,  High  Point,  North  Carolina,  The  building  consists 
of  eight  stories  and  an  unexposcd  basement;  there  is  a  mezzanine  between 
the  first  and  second  stories.  The  overall  height  of  the  building  is  about 
110  ft  and  plan  dimensions  of  50  ft  by  115  ft  provide  an  area  of  5,750 
sq  ft  on  each  floor  level.  Figure  18  shows  the  exterior  walls  and  general 
window  layout  of  the  bank.  Note  that  many  of  the  windows  on  the  first 
story  of  sides  B  and  C  have  been  bricked  In. 

The  building  has  a  structural  steel  frame  with  riveted  and  bolted 
column  and  beam  connections.  The  ribbed  floor  system  has  a  4-ln.  thick 
concrete  slab  and  4-  or  6- In.  thick  clay  tile  fillers. 

The  exterior  walls  on  sides  A  and  B  of  the  first  story  are  17-ln. 
thick  and  are  constructed  with  a  granite  veneer  and  a  brick  backing.  On 
sides  C  and  D  of  the  first  story,  the  walls  are  generally  17-ln.  thick 
solid  brick.  On  the  upper  stories,  the  walls  are  constructed  with  a 
4-ln.  thick  brick  veneer  and  an  8- In.  thick  terra  cotta  backing.  As  can 
be  noted  in  Figure  18,  the  exterior  column  lines  on  the  upper  stories  of 
sides  A  and  B  are  faced  with  a  granite  veneer.  For  all  exterior  walls, 
the  facing  Is  continuous  over  the  frame  members  and  the  backing  Is  Inset 
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COMMERCE  STREET 


SIDES  A  AND  B 


SIDES  B  AND  C 


SOUTH  MAIN  STREET 
PLAN  VIEW 

SOURCE:  RTI, 


SIDE  D 


FIGURE  18  PHOTOGRAPHS  AND  PLOT  PLAN  OF 
NORTH  CAROLINA  NATIONAL  BANK 
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In  the  frame.  The  Interior  partitions  on  the  first  story  and  mezzanine 
are  constructed  with  unreinforced  terra  cotta,  either  3-  or  6- in.  thick. 

On  the  upper  stories  the  interior  partitions  are  mostly  3-in.  unreinforced 
terra  cotta.  The  partitions  arc  nonload  bearing  and  have  numerous  open¬ 
ings  that  have  been  filled-in  with  light  wood  paneling. 

Wall  Analysis 

The  exterior  walls  on  all  sides  were  analyzed  as  unreinforced  masonry 
unit  walls  with  either  one-  or  two-way  arching.  For  Side  A  of  the  first 
story  it  was  assumed  that,  because  of  the  many  openings,  only  one-way 
arching  could  develop  between  floor  beams  on  the  first  and  mezzanine 
stories.  On  Side  B  it  was  assumed  that  one-way  arching  would  develop  in 
the  walls  between  windows.  Furthermore,  it  was  assumed  that  the  brlcked- 
In  windows  would  not  contribute  to  the  arching  strength  of  the  walls  but 
would  remain  in  place  for  a  sufficient  length  of  time  to  Influence  the 
blast  loading  and  room  filling. 

The  specific  walls  analyzed  were  as  follows: 

*  VPl:  Side  A,  wall  on  first  story;  one-way  arching  wall. 

*  VP2:  Side  B,  wall  on  first  story;  one-way  arching  wall. 

*  VP3:  All  sides,  walls  on  upper  stories;  two-way  arching  wall. 

The  results  of  the  dynamic  analysis  of  the  exterior  walls  of  the 
Bank  were: 


Predicted  Collapse  Overpressure,  psi 


Case 

Mean 

Standard 

Deviation 

10  Percent 
Probability 
Value 

90  Percent 
Probability 
Value 

VPl 

16.4 

4.2 

11.0 

21.8 

VP2 

5.4 

0.7 

4.6 

6.3 

VP3 

15.7 

4.0 

10.5 

20.8 
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Since  the  analysis  indicated  that  the  incipient  collapse  overpressure  of 
the  upper  story  walls  was  15.7  psi,  the  building  was  considered  as  a 
strong-walled  building,  and  was  therefore  a  good  candidate  for  examining 
the  dynamic  response  of  the  structural  frame. 

Frame  Analysis 

Discussion 

The  blast  loading  on  the  building  was  calculated  for  a  box-type 
building  with  openings  (Type  2,  Ref.  8),  and  the  collapse  overpressure 
of  the  exterior  walls  was  assumed  to  be  greater  than  that  of  the  frame; 
i.e.,  to  determine  the  frame  loading,  the  exterior  walls  were  assumed 
as  nonfailing.  The  blast  wave  was  assumed  to  strike  the  building  at 
normal  incidence  to  Side  B,  the  Commerce  Street  side.  The  computer 
program  determines  the  air  blast  loading  on  the  basis  of  a  specified 
peak  incident  overpressure  level  between  2  and  30  psi  for  a  10-Mt-yield 
nuclear  weapon.  Although  the  exterior  walls  had  been  analyzed  for  an 
air  blast  loading  from  a  1-Mt-yleld  weapon,  there  would  be  negligible 
difference  in  the  predicted  collapse  overpressure  of  the  walls  for  a 
10-Mt  yield. 

Three  different  types  of  frame  analyses  were  performed.  The  first 
was  an  elastic  analysis  to  determine  the  magnitude  of  the  base  shear  as 
a  function  of  the  blast  overpressure  level.  For  this  analysis,  the 
inset  exterior  walls  on  end  frame  lines  A  and  H  in  Figure  19  were  assumed 
to  act  as  shear  walls. 

The  second  and  third  types  of  frame  analyses  were  performed  to  deter¬ 
mine  the  elastic  and  inelastic  response,  respectively,  of  the  frame  acting 
alone;  for  these  analyses  the  effect  of  the  shear  wall  action  of  the  ex¬ 
terior  walls  on  the  frame  response  was  assumed  to  be  negligible.  Such 
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an  analysis  should  approximate  the  behavior  of  the  building  after  failure 
of  the  shear  walls  (l.e.,  Inset  exterior  end  walls). 

The  elastic  frame  analysis  was  conducted  at  various  Incident  over¬ 
pressure  levels  and  provided  Information  on  the  possible  collapse  strength 
of  the  frame.  Also,  the  elastic  analysis  assisted  In  the  selection  of 
the  overpressure  levels  for  running  the  more  complex  Inelastic  frame 
program. 

The  following  input  data,  required  for  the  various  analyses,  were 
obtained  directly  from  the  architectural  and  structural  plans: 

•  Percentage  of  openings  (front  face) 

•  Clearing  distance  (front  face) 

•  Floor  weight 

•  Beam  and  column  properties 

-  Moment  of  inertia 

-  Cross-sectional  area 

-  Plastic  moment  capacity 

•  Shear  wall  properties 

-  Moment  of  inertia 

-  Cross-sectional  area 

-  Plastic  moment  capacity 

-  Shear  area 

-  Shear  capacity 

Figure  19  shows  a  simplified  framing  plan  and  Figure  20  shows  the 
beam  and  column  sizes  for  frame  line  A. 
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FIGURE  19  TYPICAL  FRAMING  PLAN  (SIMPLIFIED) 
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FIGURE  20  BEAM  AND  COLUMN  SIZES  FOR 
FRAME  LINE  A 


67 


Results 


The  results  of  the  following  frame  analyses  are  summarized  in 
Tables  1  and  2: 

•  Elastic  shear  wall,  2  psi  overpressure 

•  Elastic  frame 

-  16  psi  overpressure 

-  5  psi  overpressure 

-  4  psi  overpressure 

-  3  psi  overpressure 

•  Inelastic  frame 

-  5  psi  overpressure 

-  4  psi  overpressure 

-  3  psi  overpressure 

The  maximum  stress  ratio  for  the  beams  and  columns,  shown  in  Table  1, 
is  the  ratio  of  computed  moment/yield  moment  and  for  the  shear  walls  is 
the  computed  shear/ultimate  shear  capacity.  The  ductility  ratio  listed 
in  Table  2  is  the  maximum  rotation/ yield  rotation  for  the  member.  It 
should  be  noted,  so  as  to  avoid  a  condition  of  singularity  in  the  plastic 
matrix  for  the  inelastic  frame  program,  in  Ref.  8,  "...it  is  assumed  that 
each  member  in  the  structure  is  composed  of  two  parallel  elements — one 
remains  linearly  elastic  under  all  loading  conditions  and  one  remains 
ideally  elastoplastic. . . .  For  the  computer  program  presented  in  this 
report  it  is  assumed  that  linearly  elastic  elements  represent  five  per¬ 
cent  of  the  initial  stiffness  of  the  structure.  Also,  this  makes  it 
possible  to  account  for  some  of  the  'strain  hardening'  effects  which  are 
neglected  in  the  normal  ideally  plastic  idealization." 

The  results  of  the  first  analysis,  for  the  shear  wall  building,  indi¬ 
cated  that  the  cracking  of  the  exterior  walls  on  the  ends  of  the  building 
acting  as  shear  walls  occurred  at  an  incident  overpressure  of  less  than 
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Table  1 


SUMMARY  OF  RESULTS  OF  ELASTIC  FRAME  ANALYSIS 


Maximum 

Stress  Ratio 

Axial 

Frame 

Beam 

Column 

Shear  Wall 

Force 

Line 

Moment 

Moment 

Shear 

(klps) 

Elastic.  Shear  Wall.  2  osl 

(Maximum  Deflection  - 

2.2  In.) 

A 

- 

37.3 

1.33 

0 

<1 

<1 

2.9 

COE 

<1 

<1 

6.2 

F 

<1 

<1 

9.0 

G 

<1 

<1 

5.3 

H 

- 

35,9 

2.35 

0 

Elastic.  Frame.  16  psl  (y  =  670 

—  '  . *"  'max 

1.2  In.) 

A 

43.3 

23.9 

584,4 

B 

30.8 

20.3 

886.9 

COE 

31.3 

20.2 

846.8 

F 

43.6 

20.2 

1396.6 

G 

43.6 

18.6 

968.9 

H 

33.6 

21.7 

373.6 

Elastic.  Frame.  5 

135. 

3  In.) 

A 

9.8 

4.7 

120.4 

B 

11.1 

4.0 

181.3 

COE 

6.7 

3.7 

178.9 

F 

9.8 

3.7 

292.7 

G 

9.7 

3.7 

211.2 

H 

7.8 

4.3 

77.2 

Elastic.  Frame.  4 

psl  (y 

■  max 

c  98.8 

In.) 

A 

7.2 

3.3 

87.7 

B 

8.1 

3.0 

131.9 

COE 

4.9 

3.1 

130,8 

F 

7.1 

3.1 

213.9 

G 

7.1 

2.7 

155.1 

H 

3.6 

3.2 

36.3 

Elastic.  Frame.  3 

«  63.4 

1  In.) 

A 

4.4 

2.3 

87.7 

B 

3.3 

1.6 

86,7 

COE 

3.2 

2.1 

86,4 

F 

4.6 

2.1 

141.3 

G 

4.6 

1.9 

102.8 

H 

3.7 

2.2 

37.2 

69 


Table  2 


SUMMARY  or  RESULTS  OK  INELASTIC  FRAME  ANALYSIS 


Axial 

Story 

Moment  Ratio 

Ductility  Ratio 

Force 

Level  Beam  Column  Beam  Column  (kips) 


Inelastic. 

Frame  5  psi 

(Maximum 

Deflection  =  677.6 

in. ) 

Roof 

1.91 

2.01 

27.0 

29.4 

3.7 

8 

2.51 

1.08 

29.4 

3.0 

8.9 

7 

2.51 

1.37 

30.9 

10.7 

14.1 

6 

2.69 

1.19 

35.2 

0,1 

19.7 

5 

2.72 

1.35 

35.0 

8.5 

25.3 

1 

2.63 

1.58 

32.6 

18.2 

30.7 

3 

2.31 

2.  19 

26.4 

55.7 

35.5 

2 

0.12 

1.26 

•1 

5.5 

100,6 

Mezz. 

- 

1.23 

- 

7.1 

100.6 

Inelastic , 

Frame  1  psi 

(Maximum 

Deflection  =  563.7 

in.) 

Roof 

1.81 

1.91 

23.2 

26.7 

3.5 

8 

2.25 

0.95 

25.7 

<1 

8.2 

7 

2.31 

1.32 

26.7 

9.0 

13.0 

6 

2. IS 

1.03 

29.5 

1.7 

18.1 

5 

2.45 

1.28 

29.3 

8.4 

23.2 

1 

2.31 

1.50 

26.4 

15.9 

28.1 

3 

2.03 

2.09 

21.1 

42.2 

32.4 

2 

0.09 

1.09 

<  1 

3.0 

72.6 

Mezz. 

- 

1.17 

- 

3.1 

72.6 

Inelastic  > 

Frame  3  psi 

(Max imum 

Deflection  =  249.0 

in. ) 

Roof 

1.30 

1.31 

8.4 

10.1 

2.6 

8 

1.34 

0.98 

11.8 

<1 

6.0 

7 

1.60 

1.12 

12.6 

2.6 

9.6 

6 

1.62 

0.91 

13.4 

<1 

13.2 

5 

1.64 

1.15 

13.2 

4.1 

16.8 

1.55 

1.16 

11.5 

5.8 

20.3 

3 

1.45 

1.62 

9.4 

20.6 

23.5 

2 

0.07 

1.01 

<1 

1.1 

59.3 

Mezz. 

- 

1.03 

- 

1.9 

59.3 
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2  psl,  since  the  moment  ratios  for  the  shear  walls  are  above  35  for  the 

2- psi-overpressure  level.  Therefore,  it  was  assumed  for  the  frame  anal¬ 
yses  that  the  inset  end  walls  contributed  negligible  resistance  to  the 
building,  and  that  the  analysis  of  the  frame  acting  alone  should  ade¬ 
quately  model  the  building  behavior  under  lateral  load.’*' 

An  elastic  analysis  of  the  frame  for  16  psi,  which  approximates  the 
incipient  collapse  overpressure  of  the  exterior  walls,  indicated  a  maxi¬ 
mum  stress  ratio  of  about  24  for  the  columns  and  51  for  the  beams.  Since 
the  elastic  analysis  was  much  simpler  than  the  Inelastic  analysis,  the 
frames  were  then  analyzed  for  elastic  behavior  at  5-,  4-,  and  3-psi 
overpressures  to  obtain  an  estimate  of  the  frame  strength.  The  results 
of  the  elastic  analyses,  summarized  in  Table  1,  indicate  that  the  strength 
of  the  frames  was  probably  in  the  range  of  the  lower  overpressures  exam¬ 
ined,  and  therefore  the  inelastic  frame  analyses  were  run  at  5-,  4-,  and 

3- psl  Incident  overpressure  levels. 

The  inelastic  analyses,  summarized  in  Table  2,  indicated  maximum 
ductility  ratios  at  3  psl  of  13.4  for  the  beams  and  20.6  for  the  columns, 
and  maximum  moment  ratios  of  about  1.6  for  the  beams  and  columns.  At 
the  4-pb1  level,  the  maximum  ductility  ratios  were  29.5  in  the  beams  and 
42.2  in  the  columns,  and  the  maximum  moment  ratios  were  in  excess  of  2. 

A  simplified  hand  calculation  indicated  that  the  P-A  effect,  which  is 
not  Included  in  the  computer  program,  would  increase  some  of  the  moment 


* 

Ref.  39,  which  was  a  study  of  the  effect  of  masonry  filler  walls  on 
frame  response,  indicated  that  inset  walls  do  have  a  significant  effect 
on  the  frame  resistance,  even  after  cracking  of  the  walls.  The  avail¬ 
able  frame  program  (Ref,  8)  used  could  not  model  the  effect  of  a  cracked 
inset  wall  on  the  frame  resistance  without  considerable  modification. 
However,  since  the  moment  ratios  were  quite  high  for  the  2-psl  overpres¬ 
sure  level,  and  since  only  two  of  eight  frame  lines  had  inset  walls,  it 
was  felt  that  neglecting  the  inset  walls  did  not  significantly  affect 
the  estimate  of  frame  collapse. 
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ratios  by  over  50  percent.  As  noted  on  Figure  21,  the  calculated  lateral 
deflection  of  the  top  story,  for  the  inelastic  frame  analysis,  was  about 
21  ft  for  the  3-psi  overpressure  level,  and  47  ft  for  the  4-psi  level. 

Summary 

The  results  of  the  analyses  provide  an  estimate  of  the  collapse 
strength  of  the  structural  steel  frame  of  the  bank  building  under  blast 
loading,  even  though  the  computer  program  used  cannot  predict  frame  col¬ 
lapse.  If  it  is  assumed  that  the  frame  would  collapse  at  a  ductility 
ratio  of  about  50,*  then  the  estimated  collapse  overpressure  is  between 
3-  and  4-psl  incident  overpressure  level.  The  actual  blast  strength 
could  be  much  less,  since  the  effect  of  the  axial  column  load  (P-A  effect) 
and  frame  collapse  mechanisms,  such  as  column  buckling  and  instability, 
are  not  accounted  for  in  the  computer  program. 

It  should  be  mentioned  that  the  frame  of  the  North  Carolina  National 
Bank  building  appears  to  be  constructed  of  relatively  light  structural 
shapes  that  may  not  necessarily  be  typical  of  most  NFSS  structures.  In 
any  event,  however,  the  analysis  indicated  that  the  blast  resistance  of 
the  frame  of  the  building  was  much  less  than  (possibly  only  one- fourth) 
that  of  the  exterior  walls.  This,  of  course,  is  an  important  consideration 
in  predicting  either  the  magnitude  of  building  damage  or  the  number  of 
casualties  that  might  occur  in  buildings  subjected  to  nuclear  air  blast. 


The  selection  of  a  ductility  ratio  of  50  as  indicating  a  possible  frame 
collapse  Is  arbitrary,  but  is  felt  to  be  conservative  for  estimating  a 
frame  collapse  from  the  results  of  an  inelastic  analysis.  Other  in¬ 
vestigators,  e.g..  Ref.  40,  have  used  a  frame  ductility  ratio  of  20  to 
indicate  collapse.  However,  the  actual  ductility  ratio  used  to  estimate 
collapse  is  not  too  Important,  since,  as  noted  on  Table  2,  for  an  in¬ 
crease  in  incident  overpressure  level  from  3  to  4  psi  increases  the 
maximum  column  ductility  ratio  at  story  level  three  from  20.6  to  42.2. 
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FIGURE  21  DEFLECTION  OF  FRAME  FOR  INELASTIC  ANALYSIS 
AT  3  AND  4  PSI  BLAST  OVERPRESSURE  LEVELS 
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Appendix  A 


CORRELATION  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 
FOR  LOIKilTUDINALLY  RESTRAINED  REINFORCED  CONCRETE  SLABS 
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Appendix  A 


CORRELATION  OF  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 
FOR  LONGITUDINALLY  RESTRAINED  REINFORCED  CONCRETE  SLABS 


Introduction 


This  appendix  compares  the  analytical  predictions  obtained  using 
the  theoretical  procedures  deve^'  ped  in  the  main  body  of  the  report  to 
experimental  results  from  several  tests  of  reinforced  concrete  slabs  with 
longitudinally  restrained  edges.  The  experimental  results  are  primarily 
for  statically  loaded  slabs,  although  results  for  a  few  dynamically  loaded 
slabs  are  Included.  All  slabs  were  uniformly  loaded. 

Results  from  the  following  tests  are  Included  in  the  comparison. 

•  NCEL  (Ref.  10) — Static  and  dynamic  tests  of  72-in. 
square  slabs,  3-ln.  to  6-in.  thick  (Lj =  12  to  24), 
are  reinforcement  ratios  from  0  to  1.33  percent. 

•  WES  (Ref.  11)— Static  and  dynamic  tests  of  29-in. 
square  slabs,  0.89-in.  thick  (Lj/h,  s  32.6),  and 
reinforcement  ratios  from  0,55  to  1.02  percent. 

•  MIT  (Ref.  12)— Static  tests  of  15-in.  square  slabs, 

0.75-ln.  to  1.50-in.  thick  (L;/h,  =  10  to  20),  and 
reinforcement  ratios  from  0  to  2.0  percent. 

•  Wood  (From  Ref.  10)— Static  tests  of  68-in.  square 
slabs,  2,25-ln.  thick  (Lj/h,  =  30.2),  and  reinforce¬ 
ment  ratios  from  0  to  0.26  percent. 

•  Park  (Ref.  9)— Static  tests  of  40-in.  by  60-in.  slabs 
(L^/L|  =  1.50),  0.98-ln.  to  2.0-ln.  thick  (Lj /h,  = 

20.0  to  ■10,8),  and  reinforcement  ratios  from  0  to  1.21  per¬ 
cent.  Edge  restraints  include;  (1)  all  edges  fully  re¬ 
strained;  (2)  one  short  edge  simply  supported  on  rollers, 
with  remaining  edges  fully  restrained;  and  (3)  one  long 
edge  simply  supported  on  rollers,  with  remaining  edges 
fully  restrained. 
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•  Powell  (From  Ref.  9)— Static  tests  of  20,57-in.  by 
36.0- in.  slabs  (L^/L,  =  1.75),  1.286- in.  thick, 

(Lj/h,  =  16.0),  and  reinforcement  ratios  from  0  to 
1,53  percent. 

The  results  of  these  comparisons  are  presented  in  the  following  sub¬ 
sections. 

Statically  Loaded  Slabs 

The  theoretical  resistance  functions  outlined  in  a  previous  section 
are  compared  to  results  obtained  from  static  load  tests  on  longitudinally 
restrained  reinforced  concrete  slabs.  These  comparisons  are  divided  into 
twr<  groups;  (1)  slabs  with  all  edges  fully  restrained,  and  (2)  slabs 
with  one  edge  simply  supported  and  the  remaining  edges  fully  restrained. 

Slabs  with  All  Edges  Fully  Restrained 

Dimensions,  reinforcement  data,  and  concrete  and  steel  strengths 
for  the  slabs  with  all  edges  fully  restrained  are  summarized  in  Table  A-1. 
Predicted  resistance  and  deflection  values  for  the  stages  of  behavior 
shown  in  Figure  1  in  Section  II  are  compared  with  the  test  values  in 
Table  A-2.  Average  values  of  predicted-to-test  results  are  given  for 
each  group  of  test  slabs,  as  well  as  for  the  total  number  of  slabs.  In 
general,  the  predicted  results  compare  favorably  with  the  lest  results. 

The  ultimate  flexural  resistance  computed  from  Eq.  42  is  on  the 
average  13  percent  higher  than  the  test  results.  The  greatest  disparity 
occurs  for  the  slabs  tested  by  Wood  (Kef.  10),  Park  (Ref.  9),  and  Powell 
(Ref.  9).  One  reason  for  this  may  be  that  the  concrete  strength  for 
these  slabs  is  based  on  the  cube  strength,  with  the  corresponding  cylinder 
strength,  f  ^  assumed  to  be  RO  percent  of  the  cube  strength.  This  per¬ 
centage  may  be  on  the  high  side.  The  predicted  ultimate  flexural  resis¬ 
tance  Is  also  scon  to  be  better  for  the  squore  slabs  than  for  the 
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THKORETICAL  VERSIIS  TEST  RESL'LTS  FOR  SLABS  VITH  ALL  EDCZS  FULLY  RESTRAINED 
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Tabl«  A*2  (Qoncluded) 
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rectangular  slabs.  This  may  be  due  to  the  deflection  required  to  develop 
the  crushing  strain  of  the  concrete  along  the  hinge  lines  being  different 
for  the  X-  and  z-dlrectlon  strips.  Since  the  resistance  predicted  by 
Eq.  42  Is  based  upon  the  ultimate  deflection  for  the  z-dlrectlon  strips, 
this  results  In  a  slightly  erroneous  value  for  the  x-dlrectlon  strips, 
and  thus  for  the  predicted  ultimate  resistance.  The  predicted  ultimate 
deflection  for  the  compressive  me.ibrane  region  Is  slightly  lower  than 
the  measured  values.  This  also  contributes  to  the  higher  predicted 
strengths.  As  discussed  previously,  the  ultimate  deflection  given  by 
Eq.  28  Is  bounded  by  an  empirical  upper  limit  of  0.42  times  the  slab 
thickness. 

The  predicted  secondary  resistance  Is  on  the  average  9  percent  higher 
than  the  measured  resistance.  The  corresponding  predicted  deflection  is 
on  the  average  9  percent  lower  than  the  test  results.  As  was  the  Case 
for  the  ultimate  deflection  in  the  compressive  membrane  region,  the 
secondary  deflection  is  bounded  by  an  empirical  upper  limit,  the  limit 
for  this  deflection  taken  as  three  times  the  slab  thickness.  As  Chavn  in 
subsequent  plots  of  the  load-deflection  curves,  the  secondary  resistance 
region  is  generally  not  well  defined.  This  is  particularly  true  for  the 
deflection.  The  predicted  values  are  meant  to  be  approximate  only.  In 
general.  It  Is  felt  the  difference  in  predicted  and  test  resistance  during 
this  region  has  a  minor  effect  on  the  total  predicted  behavior  of  the 
slab. 

The  predicted  values  for  the  ultimate  tensile  membrane  resistance 
and  deflection  are  both  on  the  average  approximately  25  percent  higher 
than  the  measured  values.  These  percentages  may  be  somewhat  misleading, 
however,  as  the  test  values  do  not  correspond  to  failure  values  for  all 
the  test  slabs.  This  Is  the  case  for  the  NCEL  slabs,  where  at  the  de¬ 
flections  indicated  for  y^^,  the  slabs,  though  badly  damaged,  are  still 
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Intact.  The  deflections  for  the  slabs  tested  by  Park  similarly  do  not 
indicate  collapse.  The  values  for  the  slabs  tested  at  WES  and  MIT,  how¬ 
ever,  correspond  to  significant  rupture  of  the  steel  reinforcement  and 
partial  collapse  of  the  slabs.  The  agreement  between  the  predicted  and 
test  results  is  quite  good  for  these  slabs,  particularly  for  the  collapse 
defloctlon.  A  better  value  for  comparison  purposes  may  be  the  slope  of 
the  tensile  membrane  resistance  curve  (q^^/yf^)*  The  agreement  between 
the  predicted  and  measured  slopes  of  the  tensile  membrane  resistance  is, 
on  the  average,  quite  good. 

Variation  of  the  static  resistance  with  the  center  deflection  is 
shown  for  several  of  the  test  slabs  in  Figures  A-1  through  A-4.  Both 
predicted  and  measured  results  are  given.  A  brief  discussion  of  the 
results  for  each  group  of  tests  is  also  presented  in  the  following  para¬ 
graphs. 

NCEL.  As  can  be  seen  from  Table  A-2  and  Figure  A-1,  the  predicted 
resistance  functions  for  the  NCEL  slabs  agree  quite  well  with  the  test 
results.  This  is  particularly  true  for  the  compressive  membrane  region, 
which  is  the  primary  region  of  interest  for  this  study.  In  general,  good 
agreement  is  also  found  for  the  secondary  resistance.  However,  comparison 
of  the  maximum  tensile  membrane  resistance  values  given  in  Table  A-2 
show  the  predicted  results  to  be  considerably  higher  than  the  test  results. 
As  can  be  seen  from  Figure  A-1,  the  slope  of  the  predicted  resistance 
curve  during  this  region  is  only  slightly  less  than  the  slope  of  the 
measured  resistance.  Thus,  the  discrepancy  in  the  maximum  resistance  is 
due  primarily  to  the  difference  between  the  predicted  and  measured  values 
of  the  collapse  deflection,  y^^.  The  test  values  of  y^^  for  the  NCEL 
slabs  correspond  to  rupture  of  the  steel  reinforcement  at  the  supports. 

A  review  of  the  photographs  given  in  Ref.  10  for  the  slabs  after  the 
static  tests  shows  that  although  the  slabs  are  severely  cracked,  they  are 
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FIGURE  A-1  STATIC  RESISTANCE  DIAGRAMS  FOR  NCEL  SLABS 


ANALYTICAL  RESULTS 
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FIGURE  A-1  STATIC  RESISTANCE  DIAGRAMS  FOR  NCEL  SLABS 


FIGURE  A-2  STATIC  RESISTANCE  DIAGRAMS  FOR  WES  SLABS 
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FIGURE  A-3  STATIC  RESISTANCE  DIAGRAMS  FOR  MIT  SLABS 


TEST  RESULTS 
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EDGES  FULLY  RESTRAINED 


still  intact  and  would  still  provide  protection.  How  much  reserve 
strength,  if  any,  the  slabs  still  have  is  unknown.  Thus,  comparison  of 
the  maximum  tensile  membrane  resistance  is  felt  to  be  somewhat  question¬ 
able  for  these  slabs. 

WES.  The  predicted  and  test  results  for  the  WES  slabs  arc  presented 
in  Figure  A-2.  These  results  indicate  that  the  predicted  ultimate 
strength  in  the  compressive  membrane  region  agree  quite  well  with  the 
test  results,  although  the  corresponding  predicted  ultimate  deflection 
is  about  80  percent  of  the  test  values  on  the  average.  As  indicated  in 
Table  A-1  the  steel  reinforcement  percentage  for  the  edge  strips  of  the 
slabs  was  different  than  for  the  middle  strips.  Since  the  procedures 
developed  in  this  study  assume  a  uniform  distribution  of  steel  reinforce¬ 
ment  in  each  slab  direction  (not  necessarily  the  same  in  both  directions), 
an  average  value  of  the  reinforcement  in  the  middle  and  edge  strips  was 
used  in  the  analysis.  Agreement  between  the  predicted  and  measured 
results  for  the  secondary  and  tensile  membrane  regions  is  also  quite  good. 

MIT.  The  predicted  resistance  functions  are  compared  to  the  load- 
deflection  behavior  obtained  from  tests  of  longitudinally  restrained  slabs 
conducted  at  MIT  in  Figure  A-3.  The  predicted  resistance  compares  fairly 
well  to  the  test  results  in  the  compressive  membrane  region,  but  is  some¬ 
what  high  in  the  tensile  membrane  region.  The  predicted  deflection  at 
collapse  conqiares  very  well  with  the  test  results,  however. 

Wood.  No  resistance  functions  were  plotted  for  the  slabs  tested  by 
Wood  (From  Ref.  10).  Also,  results  are  presented  in  Table  A-2  for  the 
compressive  membrane  region  only.  For  this  region,  the  predicted  resis¬ 
tance  is  approximately  25  percent  higher  than  the  test  results.  One 
possible  reason  for  this  discrepancy  may  be  the  fac,tor  used  in  converting 
the  cube  strength  of  the  concrete  to  a  corresponding  cylinder  strength. 
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The  cylinder  strength  was  assumed  to  be  80  percent  of  the  cube  strength. 
This  assumed  percentage  may  be  somewhat  higher  than  Is  actually  the  case. 

Park.  Table  A-2  compares  predicted  resistance  functions  values  to 
the  results  presented  In  Ref.  9  for  tests  on  rectangular  slabs  (length- 
to-wldth  ratio  of  1.5),  Representative  plots  of  the  predicted  and  actual 
load-deflection  curves  are  also  given  In  Figure  A-4.  The  predicted 
results  compare  fairly  well  with  the  test  results  for  the  compressive 
membrane  region.  However,  for  the  secondary  region  the  predicted  resis¬ 
tance  Is  approximately  30  percent  higher  than  the  test  results.  The 
predicted  slope  for  tensile  membrane  region  Is  considerably  below  the 
test  results.  It  should  be  noted  that  test  results  are  available  only 
for  slab  deflections  up  to  4.0  In.,  rather  than  up  to  failure.  Compari¬ 
sons  of  predicted  and  test  results  for  the  tensile  membrane  region  are 
thus  qujstlonable  for  these  slabs. 

Powell.  Results  are  presented  In  Ref.  9  only  for  the  ultimate 
compressive  membrane  strength  of  the  rectangular  slabs  (length-to-wldth 
ratio  of  1.75)  tested  by  Powell.  The  predicted  ultimate  resistance  for 
this  region  Is  on  the  average  22  percent  higher  than  the  test  results. 

Slabs  with  One  Edge  Simply  Supported 

Park  (Ref.  9),  In  addition  to  his  tests  on  rectangular  slabs  with 
all  edges  fully  restrained,  also  conducted  tests  on  two  additional  series 
of  slabs  with  one  of  their  edges  simply  supported  on  rollers,  the  remain¬ 
ing  edges  of  both  series  being  fully  restrained.  The  series  B  slabs  were 
tested  with  one  of  the  short  edges  simply  supported,  while  the  series  C 
slabs  were  tested  with  one  of  the  long  edges  simply  supported.  The  steel 
content  of  these  series  of  slabs  was  the  same  as  the  corresponding  series 
A  slab  listed  In  Table  A-1  (l.e.,  steel  contents  of  slabs  Al,  Bl,  and  Cl 
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were  identical),  except  that  there  was  no  top  steel  at  the  simply  sup¬ 
ported  edge  and  that  the  length  of  the  other  top  steel  was  such  as  to 
suit  the  particular  slab.  The  remaining  dimensions  and  properties  were 
also  the  same  except  for  the  concrete  strengths,  which  are  listed  in 
Table  A-3. 

Predicted  resistance  and  deflection  values  for  the  various  stages 
of  behavior  are  compared  with  the  test  results  in  Table  A-3.  Variation 
of  the  static  resistance  with  the  center  deflection  is  shown  for  both 
predicted  and  test  results  of  the  slabs  in  Figure  A-5.  Due  to  restric¬ 
tions  of  the  computer  program,  the  predicted  results  for  the  compressive 
membrane  region  are  based  upon  all  edges  restrained  against  rotation 
rather  than  the  actual  case  of  one  of  the  edges  simply  supported.  The 
predicted  results  reflect  the  actual  condition  of  no  longitudinal  re¬ 
straint  in  the  direction  perpendicular  to  the  simply  supported  edge, 
however. 

Because  of  the  necessity  of  treating  all  edges  as  restrained  against 
rotation  for  the  calculation  of  tha  ultimate  flexural  resistance,  the  pre¬ 
dicted  results  would  be  expected  to  be  slightly  higher  than  the  test 
results.  This  is  found  to  be  the  case  for  the  slabs  with  one  of  the  long 
edges  simply  supported:  the  predicted  ultimate  flexural  resistance  is 
37  percent  higher  on  the  average  than  the  test  results.  For  the  slabs 
with  one  of  the  short  edges  simply  supported,  however,  the  predicted 
ultimate  flexural  resistance  is  13  percent  less  than  the  test  results. 
However,  the  analytical  predictions  neglected  to  consider  the  compressive 
membrane  forces  in  the  direction  at  right  angles  to  the  simply  supported 
edges.  Some  membrane  forces  are  bound  to  develop  in  this  direction 
because  the  fully  restrained  edges  adjacent  to  the  simply  supported  edge 
are  unable  to  increase  in  length.  This  effect  is  felt  to  be  greater  for 
the  slabs  with  one  of  the  short  edges  simply  supported  due  to  the  greater 
length  of  the  fully  restrained  long  edges. 
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EDGE  SIMPLY  SUPPORTED 


An  additional  factor  in  the  conpariaon  of  the  predicted  and  actual 
ultiaate  flexural  resistance  is  the  difference  in  the  corresponding 
ultimate  deflections.  The  predicted  ultiaate  deflection  for  the  slabs 
with  one  of  the  short  edges  siaply  supported  in  0.84  in.,  equal  to  upper 
bound  empirical  limit  of  0.42  tines  the  slab  thickness.  The  correspond* 
ing  ultimate  deflection  was  only  0.30  times  the  slab  thickness,  however. 
Since  decreases  with  increasing  y^,  the  40  percent  higher  deflection 
contributes  greatly  to  the  13  percent  lower  predicted  ultiaate  flexural 
resistance.  The  opposite  effect  exists  for  the  slabs  with  one  of  the 
long  edges  siaply  supported.  The  predicted  ultiaate  deflection,  which 
for  this  case  is  calculated  using  Eq.  29,  is  on  the  average  28  percent 
lower  than  the  experlsrantal  deflection.  This  is  a  large  contributor  to 
the  37  percent  higher  predicted  ultimate  flexural  resistance,  particularly 
for  slabs  C3  and  C4.  Due  to  the  Increasing  percentage  of  steel  rela> 
forcenent  in  the  short  direction  as  compared  to  the  long  direction,  the 
value  of  3  defining  the  location  of  the  yield  lines  is  reduced,  thus 
resulting  in  reduced  values  of  the  predicted  ultiaate  deflection.  The 
predicted  ultisMte  flexural  resistance  for  these  two  slabs  is  approxi- 
Butely  70  percent  greater  than  the  corresponding  test  results. 

Test  results  for  the  secondary  resistance  are  available  only  for 
the  slabs  with  one  short  edge  slaq)ly  supported.  For  these  slabs,  the 
predicted  secondary  resistance  is  on  the  average  40  percent  higher  than 
the  test  values.  One  reason  for  this  is  again  the  treating  of  the  siaply 
supported  edge  as  being  rotationally  restrained  for  computational  pur¬ 
poses.  The  predicted  secondary  deflection  is  on  the  average  65  percent 
higher  than  the  corresponding  experimental  deflections. 

The  slope  of  the  predicted  tensile  membrane  resistance  curve 
(fift/Yfi)  tor  the  slaba  with  one  edge  simply  supported  is  approxiaately 
two-thirds  of  the  measured  results.  This  discrepancy  again,  no  doubt, 
results  from  neglecting  any  tensile  membrane  forces  in  the  long  direction 
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perpendicular  to  the  simply  supported  edge.  The  tensile  membrane  reals' 
tance  for  slab  B1  is  less  than  the  secondary  reaistance  (flexural  resis¬ 
tance  with  no  compressive  membrane  forces).  This  is  also  true  for  the 
slabs  with  one  of  the  long  edges  simply  supported.  As  was  the  case  for 
the  slabs  tested  by  Park  with  all  edges  fully  restrained,  the  measured 
deflections  for  y,^  do  not  correspond  to  collapse  of  the  slab. 

Overall,  the  predicted  results  for  the  slabs  with  one  edge  simply 
supported  and  the  remaining  edges  fully  restrained  arc  seen  to  only 
approxiutc  the  actual  results.  This  is  due,  no  doubt,  to  the  many 
assumptions  required  for  the  analysis. 

Dynamically-Loaded  Slabs 

This  section  compares  the  predicted  dynamic  behavior  of  longitudi¬ 
nally  restrained  slabs  to  the  tisie-dcf lection  response  of  dynamically 
loadet'  slabs  tested  at  NCEL  and  WES  (Refs.  10  and  11) ,  For  all  the 
slabs,  the  edges  were  fully  restrained.  Dimensions,  reinforcement  data, 
and  steel  strengths  arc  the  same  as  the  corresponding  statically  loaded 
slabs,  except  for  the  concrete  strengths,  which  are  given  in  Table  A-4. 

A  comparison  of  the  predicted  and  measured  results  for  the  dynamic  tests 
is  also  given  in  Table  A-4.  Typical  deflect  ion- time  histories,  both 
actual  and  predicted,  arc  shown  for  several  of  the  slab  tests  in 
Figures  A-6  to  A-9.  A  discussion  of  the  results  for  the  two  series  of 
testa  is  presented  in  the  following  paragraphs. 

NCEL.  Three  slabs  were  subjected  to  dynamic  loads.  Except  for  the 
concrete  strengths,  which  are  listed  in  Table  A-4,  the  slabs  were  the 
same  as  the  corresponding  statically  loaded  slabs  listed  in  Table  A-1. 

The  dynamic  loads  consisted  of  an  approximately  rectangular  load  impulse 
with  a  short  rise  time  (1  to  2  msec).  For  the  final  loading  on  each  slab, 
sudden  decreases  in  the  pressure  occurred  as  the  air  leaked  through  the 
slab. 
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Table  A-4 


THEORETICAL  VERSUS  TEST  RESULTS  FOR  DYNAMICALLY  LOAI«D  SLABS 


Peak 

Maxlaun  Dynamic 
Deflection  (ln«> 

Tine  to  Maximum 
Deflection  (nsec) 

Teat 

Load 

i££i> 

DIF 

Theory 

Theory  Test  Test 

Theory 

Theory  Test  Test 

NCEL 


3D1-1 

3.8 

3798 

40 

0.024 

0.05 

0.48 

5.2 

8.2 

1.00 

3D1-3 

10.8 

3798 

40 

0.10 

0.09 

1.11 

7.8 

7.1 

1.06 

40 

1.07 

29 

3D1-3 

39.0 

3798 

30 

5.36 

104- 

- 

55 

>100 

- 

20 

7.20 

42 

Incipient 

89.9 

3798 

40 

20 

Collapsed 

86.2 

3798 

30 

10.8-f 

m 

- 

28 

Load 

80.6 

3798 

20 

25 

4.7801 

10.5 

3320 

- 

- 

0.03 

m 

- 

4.0 

m 

4.75D3 

18.0 

3330 

40 

0.05 

0.05 

1.00 

4.0 

4.2 

0.95 

4.7801-3 

80.0 

3320 

40 

0.26 

0.37 

0.70 

6.0 

7.0 

0.86 

4.75D1-4 

62.0 

3320 

40 

0.38 

0.46 

0.83 

7.0 

7.5 

0.93 

4.78D1-5 

77.0 

3320 

40 

- 

0.60 

- 

- 

7.0 

- 

4.75D1-6 

87.0 

3320 

40 

0.76 

0.76 

1.00 

8.0 

9.4 

0,85 

40 

1.29 

30 

4.75D1-7 

98.0 

3320 

30 

1.78 

124^ 

13 

>30 

- 

20 

9.64 

10 

Incipient 

124.4 

3320 

40 

25 

Collapse 

117.0 

3320 

30 

10.84- 

- 

- 

25 

- 

ee 

Load 

109.3 

3320 

20 

25 

4.75D2-1 

11.0 

3508 

40 

0.016 

0.04 

0.40 

3.0 

4.2 

0.71 

4.75D3-2 

91.0 

3595 

40 

0.84 

0.77 

1.09 

8.0 

7.1 

1.13 

4.78D2-3 

55.0 

3898 

40 

0.32 

0.80 

0.64 

6.0 

6.0 

1.00 

4.75D2-4 

91.0 

3595 

40 

0.84 

0.96 

0.88 

8.0 

7.7 

1.04 

40 

1.72 

12 

4.75D2-8 

110.0 

3595 

30 

9.66 

9.24 

m 

30 

>28 

- 

20 

10.8 

25 

Incipient 

123.5 

3595 

40 

25 

Collapse 

116.3 

3598 

30 

10.84- 

- 

- 

28 

SB 

- 

Load 

109.4 

3598 

20 

25 

Table  A-4  (Concluded) 


Maximum  Dynamic  Time  to  Maximum 


Peak 

Deflection 

(in.) 

Deflection 

(msec) 

Test 

Load 

(psi) 

f; 

(psi) 

DIF 

(%) 

Theory 

Test 

Theory 

Test 

Theory 

Test 

Theory 

Test 

WES 

IDl 

15.0 

3683 

25 

0.48 

0.25 

1.92 

7 

IDIB 

34.7 

3683 

25 

4.35+ 

13+ 

- 

9 

24 

- 

ID2 

26.5 

3485 

25 

4.35+ 

7.5+ 

0.58 

12 

31 

0.39 

IDS 

27.8 

3485 

25 

4.35+ 

13+ 

- 

12 

31 

0.39 

ID4 

27.0 

3485 

25 

4.35+ 

7.0 

0.62 

12 

34 

- 

IDS 

27.4 

3485 

25 

4.35+ 

13+ 

- 

12 

15 

- 

Incipient 

18.2 

3683 

25 

4.35+ 

- 

- 

20 

- 

- 

Collapse 

Load 

26.4* 

3683 

25 

•  7.25+ 

16 

I  IDl 

56.8 

3810 

25 

4.35+ 

13+ 

- 

7 

15 

- 

IID2 

36.6 

3810 

25 

4.35+ 

8.1 

0.54 

10 

28 

0.36 

I  IDS 

43.1 

3810 

25 

4.35+ 

13+ 

- 

9 

28 

- 

Incipient 

25.6 

3810 

25 

4.35+ 

- 

- 

14 

- 

- 

Collapse 

Load 

38.8 

3810 

25 

7.25+ 

- 

- 

14 

- 

- 

*  Based  on  collapse  defl<°ction  of  0.25  of  the  slab  length. 
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FIGURE  A-8  DYNAMIC  DEFLECTION-TIME  HISTORIES  FOR  NCEL  SLAB  4.75D2 
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FIGURE  A>d  DYNAMIC  DEFLECTION-TIME  HISTORIES  FOR  WES  SLABS 
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Slab  3D1  required  three  cycles  of  loading  to  failure,  while  slabs 
4.75D1  and  4.75D2  were  subjected  to  seven  and  five  cycles  of  loading, 
respectively.  All  slabs  received  very  little  damage  prior  to  the  last 
cycle  of  loading.  The  last  cycle  of  loaaing  caused  a  local  failure  of 
slab  3Di  and  total  destruction  of  the  4.75-ln.  thick  slabs. 

Results  in  Ref.  10  indicate  a  dynamic  increase  factor  (DIF)  of  30 
to  60  percent  for  the  dynamically  loaded  slabs.  A  value  of  40  percent 
was  used  for  the  predicted  results  given  in  Table  A-4.  Except  for  the 
last  cycle  of  loading  for  each  slab,  the  predicted  maximum  deflection  and 
corresponding  time  of  occurrence  are  in  good  agreement  with  the  test 
results.  For  the  last  cycle  of  loading,  however,  the  predicted  results 
do  not  indicate  failure,  as  occurred  during  the  tests.  Therefore,  addi¬ 
tional  predictions  were  made  for  20-  and  30-percent  dynamic  increase 
factors.  Predicted  time-deflection  histories  arc  compared  with  test 
results  for  several  of  these  tests  in  Figures  A-6  to  A-8.  Results  are 
discussed  for  each  of  the  slabs  in  the  following  paragraphs. 

The  predicted  deflection  time-histories  shown  for  test  3D1-3  in 
Figure  A-6  indicate  that  the  best  agreement  with  test  results  occurs 
for  a  dynamic  increase  factor  of  30  percent.  As  Indicated  in  Table  A-4, 
however,  the  best  agreement  for  maximum  deflection  occurs  for  a  DIF  of 
20  percent.  This  discrepancy  may  be  due,  at  least  in  part,  to  the  local 
failure  experienced  by  slab  3D1. 

Predicted  deflectlon-tlme  histories  are  compared  to  the  results 
from  tests  4.75D1-2,  4.75D1-4,  and  4.7SD1-7  in  Figure  A-7.  The  predicted 
deflection-time  curve  for  test  4.7SD1-2  agrees  quite  well  with  the  test 
results.  For  test  4.75D1-4,  the  predicted  maximum  deflection  is  somewhat 
‘'css  than  that  experienced  during  the  test.  After  the  maximum  deflection 
occurs,  the  predicted  deflect  ion- time  curve  deviates  considerably  from 
the  test  results.  This  large  deviation  results  from  neglecting  damping 
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In  the  dynamic  analysis.  For  test  4.75D1-7,  the  best  agreement  between 
predicted  and  actual  results  occurs  for  a  dynamic  increase  factor  of 
20  percent. 

Predicted  deflectlon-tlme  histories  are  compared  to  results  from 
tests  4.75D2-2  and  4.7502-3  In  Figure  A-8.  The  predicted  deflections 
are  In  fairly  good  agreement  with  the  actual  results  up  to  maximum  de¬ 
flection.  After  this,  again  because  damping  was  neglected  In  the  dynamic 
analysis,  the  predicted  deflections  deviate  widely  from  the  actual  results. 
No  results  were  plotted  for  test  4.75D2-5.  The  results  from  Table  A-4, 
however,  indicate  the  best  agreement  between  predicted  and  test  results 
occurs  for  a  dynamic  Increase  factor  of  30  percent. 

Incipient  collapse  pressures  were  also  predicted  for  the  three  test 
slabs,  using  dynamic  increase  factor  factors  of  20,  30,  and  40  percent. 
These  values  are  shown  In  Table  A-4.  Comparison  of  these  values  with  the 
measured  peak  loads  for  the  last  cycle  of  loading  for  each  slab  indicates 
that  the  best  agreement  occurs  for  a  20-percent  dynamic  increase  factor. 
Even  these  values  are  slightly  higher  than  the  actual  test  results.  Pos¬ 
sible  reasons  for  this  may  be  inaccuracies  in  the  predicted  resistance 
function,  particularly  for  the  tensile  membrane  region,  and  possible 
damage  experienced  by  the  slabs  during  the  prior  test  load  cycles. 

WES.  Dynamic  tests  were  conducted  on  eight  slabs.  Properties  for 
these  slabs  were  the  same  as  listed  for  the  corresponding  series  I  and 
II  slabs  In  Table  A-1,  except  for  the  concrete  strengths,  which  are 
listed  in  Table  A-4.  A  comparison  of  the  predicted  and  test  results  is 
also  given  in  Table  A-4.  The  dynamic  loads  consisted  of  an  approximate 
triangular  Impulse  with  a  rise  time  of  about  4  msec,  followed  by  a  fairly 
rapid  decay  until  about  20  msec,  and  a  slower  decrease  to  zero  at  about 
100  to  120  msec.  Except  for  test  IDl,  the  pressures  were  such  that  col¬ 
lapse  occurred  on  the  first  cycle  of  loading.  Due  to  bad  primacord,  a 
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low  pressure  of  15  psi  resulted  for  test  IDl,  with  only  hairline  cracks  ^ 
occurring  around  the  periphery  of  the  slab.  A  second  test  load,  IDIB, 
collapsed  the  slabs. 

Based  on  results  from  Ref.  11,  a  dynamic  increase  factor  of  25  per¬ 
cent  was  used  for  the  predicted  results  given  in  Table  A-4.  A  comparison 
of  these  results  with  the  test  results  at  first  glance  indicates  a  wide 
disparity.  It  should  be  noted,  however,  that  the  values  given  for  the 
maximum  test  deflections  include  behavior  after  collapse  of  the  slab  has 
occurred,  while  the  predicted  results  are  only  concerned  with  behavior 
up  to  the  predicted  collapse  deflection.  Since  so  many  of  the  test  slabs 
experienced  total  collapse,  comparisons  between  the  predicted  and  test 
results  a. e  inconclusive.  However,  disregarding  deflections  of  the  slabs 
that  totally  collapsed  (and  also  test  IDl),  the  maximum  deflections  aver¬ 
aged  7.25  in.  (25  percent  of  the  slab  span)  for  the  series  I  slabs  and 
8.1  in.  (28  percent  of  the  slab  span)  for  the  series  II  slabs.  The  pre¬ 
dicted  results,  however,  indicate  collapse  to  occur  at  4.35  in.,  or  50 
to  60  percent  of  the  actual  deflections.  The  reason  for  this  discrepancy 
may  be  due  to  the  fact  that  statically-loaded  slabs  (on  which  the  pre¬ 
dicted  resistance  is  based)  start  to  fall  in  one  area,  subsequently 
unloading  the  rest  of  the  slab.  Thus,  failure  is  concentrated  in  one 
area  of  a  statically  loaded  slab.  Dynamically  loaded  slabs,  however, 
tend  to  fall  more  uniformly,  with  the  quicker  response  of  the  slab 
evening  out  the  local  defects.  This  may  result  in  a  larger  deflection 
since  local  effects  are  minimized. 

Predicted  and  actual  deflection-time  histories  are  shown  for  several 
of  the  slabs  in  Figure  A-9.  These  results  indicate  that,  for  the  periods 
of  comparison,  the  predicted  response  is  quicker  than  actually  occurred. 
One  reason  for  this  difference  is  the  use  of  the  load-mass  factor  (used 
in  simplifying  the  actual  slab  to  an  equivalent  single  degree-of-freedom 
system)  corresponding  to  the  plastic  region  for  the  entire  range  of 
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behavior.  During  the  Initial  stages  of  behavior,  the  actual  slab  Is 
somewhat  stlffer  than  Indicated  by  use  of  this  factor.  While  having  a 
minor  effect  on  the  maximum  response,  the  effect  on  the  time  to  maximum 
response  Is  less  than  actually  occurred.  This  effect  Is  also  observed 
In  comparing  the  predicted  and  actual  time  to  maximum  deflection  In 
Table  A-4.  This  effect  Is  felt  to  be  minor  as  far  as  the  Incipient 
collapse  overpressure  Is  concerned. 

Predicted  Incipient  collapse  overpressures  for  the  series  I  and  II 
slabs  were  calculated  as  18.2  psl  and  25.6  psl,  respectively.  These 
values  correspond  to  the  upper  bound  average  values  of  28.7  psl  and  45.5 
psl,  respectively,  obtained  from  the  tests.  Modifying  the  resistance 
functions  to  reflect  the  collapse  deflections  observed  In  the  test  (0.25 
of  the  slab  length)  Increases  the  predicted  Incipient  collapse  overpres¬ 
sures  to  26.4  psl  and  38.8  psl  for  the  series  I  and  II  slabs,  respectively. 
These  correspond  quite  favorably  to  the  actual  results. 

Summary 

In  general,  the  analytical  predictions  for  the  static  resistance  of 
longitudinally  restrained,  reinforced  concrete  slabs  were  found  to  compare 
favorably  with  experimental  results  from  several  tests  of  statically  and 
dynamically  loaded  slabs.  This  correlation  was  extremely  good  for  the 
compressive  membrane  region  for  the  square  slabs  fully  restrained  on  all 
four  edges.  The  predicted  ultimate  compressive  membrane  resistance  for 
rectangular  slabs  fully  restrained  on  all  edges  was  found  to  be  slightly 
on  the  high  side.  Predicted  secondary  resistance  and  deflection  values 
were  also  found  to  compare  favorably  with  experimental  results.  Although 
predicted  values  for  both  the  ultimate  tensile  membrane  resistance  and 
deflection  were  about  25  percent  higher  than  the  measured  results,  these 
comparisons  are  somewhat  misleading,  since  not  all  of  the  slabs  were 
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tested  to  collapse.  Agreement  between  predicted  and  measured  results 
for  the  slope  of  the  tensile  membrane  resistance  was  quite  good. 

A  limited  number  of  comparisons  were  also  made  for  rectangular 
slabs  with  one  edge  simply  supported  on  rollers  and  the  remaining  edges 
fully  restrained.  In  general,  correlation  between  the  predicted  and 
measured  resistance  function  values  was  only  fair.  These  results  would 
tend  to  Indicate  that  for  such  slabs  the  predicted  resistance  function 
should  be  considered  as  approximate  only. 

Predicted  dynamic  response  was  also  compared  with  the  time-deflect Ion 
behavior  of  several  dynamically  loaded  square  slabs  with  all  edges  longi¬ 
tudinally  restrained.  The  slab  response  model  generally  provided  an 
adequate  prediction  of  the  experimental  behavior.  Dynamic  Increase  fac¬ 
tors  of  20  to  40  percent  were  used  In  calculating  the  predicted  response. 
Best  correlations  were  obtained  for  a  dynamic  Increase  factor  of  about 
25  percent.  Experimental  results  also  Indicated  d)mamlc  collapse  deflec¬ 
tions  on  the  order  of  25  percent  of  the  slab  length,  as  compared  to  values 
of  approximately  15  percent  of  the  slab  length  obtained  from  the  static 
tests  used  In  most  of  the  analyses.  This  Increase  Is  felt  to  be  due  to 
tho  minimizing  of  the  effect  of  local  failures  as  a  result  of  the  rapid 
response  of  the  dynamically  loaded  slabs  as  compared  to  the  statically 
loaded  slabs. 
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NmtENCLATURE 


A  Area  of  section 

A,  Area  of  tension  steel  In  reinforced  concrete 

slab  per  unit  width,  sq  In. /In. 

A,'  Area  of  compression  steel  In  reinforced 

concrete  slab  per  unit  width,  sq  In. /in. 

A«l  Area  of  reinforcing  steel  per  unit  width 

continuous  along  the  long  span  In  reinforced 
concrete  slab,  sq  In. /In. 

At  s  Area  of  reinforcing  steel  per  unit  width 

continuous  along  short  span  In  reinforced 
concrete  slab,  sq  In. '’In. 

b  Width  of  cross  section.  In. 

Cf  Coefficient  relating  effective  nnment  of 

Inertia  of  wood-Jolst  floor  and  nrnnent  of 
Inertia  of  joist 

d  Distance  from  compressive  face  of  reinforced 

concrete  slab  to  centroid  of  tension  steel.  In. 

d'  Distance  from  compressive  face  of  reinforced 

concrete  slnb  to  centroid  of  compressive 
steel.  In. 

E  Modulus  of  elasticity,  psl 

E]  Modulus  of  elasticity  of  wood  Joist,  psl 


f,  Stress  In  tension  steel, 

t'  Stress  In  compression  ate 

f.  Unit  stress  In  wood,  psl 

f„,,  Proportional  limit  In  com 

f«eu  Ultimate  unit  compressive 

fvt  Unit  tensile  stress  In  wo 

fktu  Ultimate  unit  tensile  str 

fy  Static  yield  strength  of 

hj  Height  of  wood  Joist,  In. 

h|  Thickness  of  concrete  sla 

I  Moment  of  Inertia,  In.* 

I,  Effective  moment  of  Inert 

floor.  In.* 

Ij  Moment  of  Inertia  of  wood 

kt ,kt  Tensile  membrane  factors 

ku  Coefficient  Indicating  po 

axis  of  reinforced  concre 
strength 

kub  Coefficient  Indicating  po 

axis  of  reinforced  concre 
taneous  yielding  of  tens! 
of  concrete  (balanced  con 


E,  Moduli s  of  elasticity  of  steel,  psl 

Et  Total  energy  Input  of  uniform  load  acting  on 

slab.  In. -lb 

E,  ,Et  Energy  Input  of  uniform  load  for  slab  strips 
In  tho  x-dlrectlon  and  z-dlrectlon,  lb 

(El)'  Quantity  defined  by  Eqs.  70  and  73 

(El).ri  Quantity  defined  by  Eq.  71 

fe  Compressive  strength  of  6-  by  12-ln.  concrete 

cylinder,  psl 

fe"  Compressive  strength  of  concrete  In  flexure, 

(kafp),  psl 

tic  Dynamic  compressive  strength  of  concrete,  psl 


ki  Coefficient  relating  the 

compressive  stress  to  the 
strength  of  concrete  In  f 

ka  Coefficient  relating  the 

strength  of  concrete  In  f 
presslve  strength  of  6-  b 

Klm  Load-mass  factor 

L  Span  length.  In. 

Li  Length  of  slab  In  long  dl 

Ls  Length  of  slab  In  short  d 

M  Bending  moment.  In. -lb 

Me  Bending  moment  at  center 


tiy  Dynamic  yield  strength  of  reinforcing  steel, 

psl 

f.  Modulus  of  rupture,  psl 


Mu  Ultimate  moment  capacity 

mldthlckness  of  slab.  In. 
Ultimate  moment  capacity 
Joist,  In. -lb 


Pracediii  mi  Milk 
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f.  Stress  In  tension  steel,  psi 

t"  Stress  In  compression  steel,  psi 

Unit  stress  in  wood,  psi 

Proportional  limit  In  compression  In  wood,  psi 

Ultimate  unit  compressive  stress  in  wood,  psi 

Unit  tensile  stress  in  wood,  psi 

f«tu  Ultimate  unit  tensile  stress  in  wood,  psi 

fy  Static  yield  strength  of  reinforcing  steel,  psi 

hj  Height  of  wood  Joist,  in. 

h.  Thickness  of  concrete  slab,  in. 

I  Moment  of  Inertia,  in.* 

I,  Effective  moment  of  inertia  of  wood-joist 

floor,  in.* 

Ij  Moment  of  inertia  of  wood  Joist,  in.* 

kt ,k(  Tensile  membrane  factors 

ku  Coefficient  indicating  position  of  neutral 

axis  of  reinforced  concrete  slab  at  ultimate 
strength 

kub  Coefficient  indicating  position  of  neutral 

axis  of  reinforced  concrete  slab  at  simul¬ 
taneous  yielding  of  tension  steel  and  crushing 
of  concrete  (balanced  condition) 

ki  Coefficient  relating  the  average  concrete 

compressive  stress  to  the  maximum  compressive 
strength  of  concrete  in  flexure 

ks  Coefficient  relating  the  maximum  compressive 

strength  of  concrete  in  flexure  to  the  com¬ 
pressive  strength  of  6-  by  12-in.  cylinders 

Kim  Ltad-mass  factor 

L  Span  length,  in. 

L|  Length  of  slab  in  long  direction,  in. 

Ls  Length  of  slab  in  short  direction,  in. 

M  Bending  moment,  in. -lb 

He  Bending  moment  at  center  of  span,  in. -lb 

Mu  Ultimate  moment  capacity  per  unit  width  about 

mldthlcknesB  of  slab,  in.-lb/in. 

Ultimate  moment  capacity  at  center  of  wood 
Joist,  in. -lb 
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Ultimate  moment  capacity  per  unit  width  about 
mid  thickness  of  section  under  the  average 
thru.it  along  the  span,  In.-lb/ln. 

Thrui^t  per  unit  width  causing  simultaneous 
yielding  of  tension  s'-eel  and  crushing  of 
concrete  at  section,  lb/ In. 

Thrust  per  unit  width  at  ultimate  behavior 
(Initial  crushing  of  extreme  fiber  of  concrete 
at  hinge  section),  Ib/ln. 

,Nut  Thrust  per  unit  width  at  ultimate  behavior 

for  slab  strips  In  x-dlrectlon  and  z-dlrectlon, 
Ib/ln. 

Ratio  of  ultimate  tensile  stress  and  ultimate 
compressive  stress  of  wood  beam 

Steel  ratio,  tension  steel 

Steel  ratio,  compression  steel 

Peak  incident  overpressure,  psl 

Unit  resistance  for  uniformly  loaded  member, 
psl 

Tensile  membrane  resistance  per  unit  area  of 
slab  corresponding  to  yft,  psl 

Secondary  resistance  per  unit  width  of  longi¬ 
tudinally  restrained  reinforced  concrete  slab, 
pal 

Ultimate  flexural  resistance  per  unit  area  of 
longitudinally  restralnded  reinforced  concrete 
slab,  psl 

Ultimate  resistance  per  unit  area  of  wood- 
Jolst  floor  system,  psi 

Maximum  resistance  per  unit  area  for  elastic 
phase,  psl 

Reinforcing  Index,  d(p  -  p')f,/fe",  in. 
Reinforcing  index,  d(pf,  -  p'f')/f",  in. 

Joist  spacing,  In. 

Ratio  of  longitudinal  movement  of  slab  in  x- 
dlrectlon  to  slab  length  In  long  direction 
resulting  from  partial  restraint 

Ratio  of  longitudinal  movement  of  slab  in  z- 
directlon  to  slab  length  In  short  direction 
resulting  from  partial  restraint 


•••■i 


Wt 


X 

y 

yn 


yi 


yt 


y» 


yi 

z 

or 


a 

Yx.y* 

c 

C. 

*• 

«u 

e, 

p 

9 


Longitudinal  movement  per  unit  la 
strip  in  x-dlrectlon  and  z-dlrec1^ 
from  effect  of  partial  restraint  j 
strain,  in. /in. 

’i  >tal  work  done  by  the  resisting  t 
thrust  on  the  slab,  in. -lb 

Work  done  by  the  resisting  moment 
on  the  slab  strips  in  the  x-dlreoi 
direction,  lb 

Distance  along  slab  in  long  dlreej 

Deflection  of  slab  along  yield  11) 

Deflection  at  center  of  beam  or  s] 

Deflection  at  collapse  of  slab  in 
membrane  mode,  in. 

Deflection  at  development  of  secoi 
tance  of  slab,  in. 

Deflection  at  intersection  of  seen 
tensile  membrane  phases,  in. 

Deflection  at  ultimate  resistance 
beam,  in. 

Maximum  deflection  for  elastic  phi 

Distance  along  slab  in  short  diret 

Coefficient  relating  the  height  ol 
to  the  distance  from  the  compress! 
the  neutral  axis 

Coefficient  defining  location  of  j 
Parameters  defined  by  Eqs.  la  and 
Strain,  in. /in. 

Total  axial  strain  along  a  slab  st 
from  the  compressive  membrane  fore 

Strain  in  tension  steel,  in. /in. 

Strain  in  compression  steel,  in./l 

Ultimate  concrete  strain  in  flexui 

Angle  of  rotation  of  slab  strip  in 
radians 

Radius  of  curvature  of  neutral  axl 
beam 

Slope  of  end  portion  of  slab  strip 


Pnnffii  mil  unk 
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I|>«city  per  unit  width  about 
let  ion  under  the  average 
ipan,  in.-lb/in. 

dth  cauaing  simultaneous 
m  steel  and  crushing  of 
■  t  lb  in. 

irtth  at  ultimate  behavior 
of  extreme  fiber  of  concrete 
Ib/in. 

«th  at  ultimate  behavior 
;  x-directlon  and  s-direction, 


tensile  stress  and  ultimate 
.  of  wood  beam 

steel 

^slon  steel 
ressuro,  psl 

r  uniformly  loaded  member, 

Mlstance  per  urit  area  of 
to  Xft  .  psi 

»  per  unit  width  of  longl- 
I  reinforced  concrete  slab, 

Ml stance  per  unit  area  of 
painded  reinforced  concrete 

per  unit  area  of  wood- 

»el 

•r  unit  area  for  elastic 

IP  -  p'lVC,  in. 

P'.  -  P'f;)/C.  in. 


movement  of  slab  in  x- 
gth  in  long  direction 
I  restraint 

movement  of  slab  in  z- 
fth  in  short  direction 
I  restraint 


Longitudinal  movement  per  unit  length  of  slab 
strip  in  x-dlrectlon  and  z-directlon  resulting 
from  effect  of  partial  restraint  and  axial 
strain,  in. /in. 

Wt  Total  work  done  by  the  resisting  moment  and 

thrust  on  the  slab,  in.^lb 

Vi,W,  Work  done  by  the  resisting  moment  and  thrust 
on  the  slab  strips  In  the  x-directlon  and  z- 
direction,  lb 

X  Distance  along  slab  in  long  direction,  In. 

y  Deflection  of  slab  along  yield  line.  In. 

y,  Deflection  at  center  of  beam  or  slab,  in. 

fft  Deflection  at  collapse  of  slab  in  tensile 

membrane  mode,  in. 

y,  Deflection  at  development  of  secondary  resis¬ 

tance  of  slab,  in. 

yx  Deflection  at  intersection  of  secondary  and 

tensile  membrane  phases,  in. 

yv  Deflection  at  ultimate  resistance  of  slab  or 

beam,  in. 

yi  Maximum  deflection  for  elastic  phase,  in. 

z  Distance  along  slab  In  short  direction,  in. 

a  Coefficient  relating  the  height  of  a  wood  beam 

to  the  distance  from  the  compression  face  to 
the  neutral  axis 

B  Coefficient  defining  location  of  yield  lines 

yi,ya  Parameters  defined  by  Eqs.  la  and  lb 
c  Strain,  in. /in. 

C(  Total  axial  strain  along  a  slab  strip  resulting 

from  the  compressive  membrane  force,  in. /in. 

c,  Strain  in  tension  steel,  in. /in. 

Strain  in  compression  steel,  in. /in. 

(u  Ultimate  concrete  strain  in  flexure,  in. /in. 

Ox  Angle  of  rotation  of  slab  strip  in  x-dlrection, 

rad ians 

p  Radius  of  curvature  of  neutral  axis  of  bending 

beam 

9  Slope  of  end  portion  of  slab  strip,  radians 
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